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ABSTRACT 



In this thesis several incompressible oscillatory flow and 
flutter problems were investigated. First, a previously 
developed unsteady panel code was modified so that systematic 
comparisons with Theodorsen's classical theory could be 
accomplished. It was found that the panel code is in excellent 
agreement with the Theodorsen results. Second, the panel code 
was applied to the analysis of bending- torsion flutter. Again, 
general agreement with Theodorsen's flutter predictions was 
obtained. In the experimental part of the thesis two flow 
visualization experiments were performed. First, the vortical 
flow patterns generated by an airfoil executing harmonic 
plunge oscillations were visualized. In the second experiment, 
the interference effects between a stationary airfoil and a 
small vane executing plunge oscillations were explored. 
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I . INTRODUCTION 



A. GENERAL 

In this thesis, several numerical methods were used to 
analyze the flow about an airfoil performing unsteady motion 
in an inviscid incompressible fluid. First, the unsteady 
motion of a single airfoil was analyzed after modifying the 
U2DIIF code [ref. 2]. The primary purpose was to verify the 
code against the proven theory of Theodorsen for analyzing the 
phenomenon of flutter. To accomplish this the U2DIIF code was 
modified to calculate aerodynamic values over a range of 
reduced frequencies and then apply these values to the 
flutter analysis. 

Next, the propulsive effects of a plunging airfoil were 
verified through experimental methods using a low speed 
plexiglas wind tunnel. 

Finally, an exploratory test was conducted in the 
department's smoke tunnel to study the interaction between a 
plunging airfoil and a stationary large airfoil. 

B . SCOPE 

Chapter II contains the modification of the single airfoil 
U2DIIF code into the code UPOT.f and extensive verification of 
this code against results produced by Theodorsen. Chapter III 
describes the UPOT code and explains the modifications which 
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were added to solve the flutter determinant. In chapter IV the 
flow visualization experiment is described which was performed 
to study the vortical wake patterns produced by a plunging 
airfoil. In chapter V a second experiment is described which 
was performed to explore a plunging airfoil's potential for 
control of flow separation. 
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II. SINGLE AIRFOIL ANALYSIS 



A. U2DIFF PANEL CODE 

1 . Geometry 

Figure 2.1 shows a representation of the system that 
is analyzed using the panel code. Shown are the values for h 
(plunge) and oi (AOA) . 

2. U2DIFF 

The U2DIIF code was developed by TENG [ref .2] for the 
study of unsteady inviscid and incompressible flow over a 
single airfoil. The code is based on the extension of the 
panel method, developed by Hess & Smith [ref. 4] for steady 
potential flow problems, to include the unsteady motion of the 
airfoil that is continuously shedding vortices into the 
trailing wake. This vortex shedding process is nonlinear in 
that the wake vortices influence the flow over the airfoil 
which in turn alters the vortex shedding as the airfoil 
proceeds in time. 

The non-linearity of the unsteady flow makes this 
problem different from the steady flow problem which requires 
only simple Gaussian elimination. Teng developed a code that 
used an iterative type of solution. Typical program output 
includes the airfoil pressure distribution, force and moment 
coefficients, and the trailing vortex wake pattern. No 
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Figure 2.1 Airfoil Geometry 



4 



attempt is made here to reproduce the work of Teng or to 
explore the operation of the U2DIIF code, but the reader is 
encouraged to review reference 2 . 

B . PHASE PROGRAM 

The phase program was put through some verification by 
Neace [ref. 9] and modified slightly in order to present 
results for harmonic motion. The code PHV3 . f (phaseshift) was 
written by Neace to convert the time dependent output of lift 
and moment histories to harmonic output using an iterative 
curve fit algorithm: 



Fit) =Amp*Sin{(jit+^) (2.1) 

where Amp = amplitude of motion, u = frequency, and 0 = phase 
angle between motion and the aerodynamic forces. One primary 
output of this program was the values of phaseshift (0) 
between the AOA and coefficients of lift (C^) and moment (C„) 
for the pitching airfoil and the phaseshift between the 
plunge value (h/2b) and the C,^ and C„ for the plunging airfoil . 
The other output was the amplitude of and C„ for the 
pitching or plunging case. 

C. MODIFICATION OF U2DIIF AND PHASE PROGRAM 

In an attempt to make the above mentioned codes more "user 
friendly", the two codes were combined into a single code 
named UPOT.f. The modification involved a new input file 
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called UPOT.in which gives the user several options of 
operation. The input file can call for the analysis of steady 
flow only, straight and modified ramp motion, pitch 
oscillation, plunge oscillation, and the capability of 
performing the oscillation analysis over a series of reduced 
frequencies. A sample input file is shown in Figure (2.2). 

1 . Output 

Outputs from the code have been limited to reduce the 
amount of computer space taken up by the code operation. A 
sample calculation was run using the input from Figure (2.2) 
and the on-screen output is shown in Figure (2.3) . Of course, 
the user can modify the output portions of the code to 
minimize output. The following list describes the input/output 
files and the data they contain. 

a. UPOT.IN; The input file figure (2.2). 

b. CL.d: This file contains the various AOA values along 

with its corresponding for each time step. 

c. CM.D: This file contains the various AOA values along 

with its corresponding C„ for each time step. 

d. PHASE, d: This file contains the values of non 

dimensional time (t) , AOA, C^, C„, for each time step. 

e. FOR015.DAT: This file contains the values of non 

dimensional AOA, curve fit for C^, curve fit for C^, 
(used in the phase portion of program) . 

f. CPSS.d: This file contains the steady state pressure 

coefficient for the mid point locations of all the air 
foil panels. 

g. CPUOOS.d: This file contains the unsteady pressure 

coefficient for the mid point location of all airfoil 
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panels, (in this case the values are for an AOA equal 
to 5 degrees) . 

h. PHZSWP.d: This file contains the phase information of 

the reduced frequency sweep portion for the program. 
The file contains the phase angle of C^, and C„, and 
the amplitude of C^, C„. 

i. FLUTTER. IN: This file contains information that can 

be used to solve the flutter determinant. It contains 
Kp, Re, Im, C„ Re, C„ Im for the pitch or plunge 
case . 
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AIRFOIL TYPE : NACA 0012 AIRFOIL 
NLOWER = 50 , NUPPER = 50 



I FLAG NLOWER 
0 50 

AIRFOIL TYPE 
7 



IRAMP 

0 

FREQ 

.68 

IGUST 

0 

ITRANS 

0 

CYCLE 

2 

naoc 

2 



lOSCIL 

1 

RFQSTP 

O.Cl 

UGUST 

0 . 

DELHX 

.00 

NTCYCLE 

60 

& naot X 
05 10 



NUPPER 

50 



ALP I ALPMAX PIVOT 

-3.0 3.0 0.37 

RPQFNL 

.7 

VGUST 

0 . 

DELHY DELI 
.0 .0 
TOL 
0.005 

aoa values multiplied by 1C (integer) 
20 25 39 50 



PHASE 

0.00 



IRAMP 0: n/a RFREQ is based on fol. chora 

1 : St ra iqht ramp 
2: Modified ramp 

lOSCIL 0: n/a RFREQ is basea on full chord 

1: Sinusoidal pitch, motion starts at nin Aoa 

ITRANS 0: n/a 

1: Translational harmonic oscillation 

ALPI/ALPMAX Minimum/MAX AOA in degrees for I RAMP/ ITRANS/ lOSCI L 
MAX does not apply for ITRANS 

PIVOT Location of Elastic Axis as a fraction of fuil chord 
FREQ Initial reduced frequency for program 

RFQSTP Reduced freq step size for a sweep of f req .' s (enter 0.0 if only one calcu 
lation is desired.) 

RFQFNL Final freq for the sweep 

DELHX Translational amount in the chordwise direction (dist/full chord) 

DELHY Max Translational amount in the vertical direction ( h/ f ul Ichord (b) } 

DELI Min Translational amount in the vertical direction (h/b) 

CYCLE : • of cycles for oscillatory motions 

-In case of ramp, cycle=1.5 denotes airfoil is held 
at max aoa for the duration of .5 cycle 
-For steady state solution set it to 0 

NTCYCLE: • of time steps for each cycie 

CYCLE'NTCYCLE is limited to 200 currently. 

TOL Tolerance for convergence of the unsteady solut ion .( recommend using not 

iess than .001) 

NAOT: f of input aoa for cp output 

- angles should be in increasing order, 

- for oscilatory motions angles should increase 
first, then decrease. Decreasing angles are for 
the return cycle.. 



Figure 2.2 UPOT.IN 
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stdin 


Page 1 


4 




AIRFOIL TYPE : NACA 0012 , 
NLOWER - 50 , NUPPER = 50 


AIRFCIL 












I FLAG (C:NACA, 1: INPUT) 


0 




NO. 


PANELS UPPER SURFACE 


50 




NO. 


PANELS LOWER SURFACE 


50 




7 

OSCILLATORY MOTION, lOSCIL 
INITIAL ANGLE OF ATTACK 


-3. 0000 




FINAL ANGLE OF ATTACK 


3.0000 




REDUCED FREQ. FOR OSCIL 


0.6800 




REDUCED FREQ. STEP = 0.0100 




FINAL REDUCED FREQ. 


0.7000 




PIVOT POINT 


0.3700 




TOTA 


L i OF CYCLES 


2.0000 




1 of 


TIME STEPS PER CYCLE = 


60 




TOLERANCE FOR CONVERGENCE 


0.0050 




FREQ SWEEP 






FREQ ^ 


0.680000 






STEADY FLOW SOLUTION AT ALPHA = - 


3.000000 




1 


0.999507 -1.116717 -0.086820 


0.280753 -0.848084 




2 


0.997535 -1.108171 -0.086820 


0.213074 -0.887088 




3 


0.993600 -1.100064 -0.086820 


0.168200 -0.912031 




4 


0.987718 -1.092523 -0.086820 


0.134381 -0.930386 




5 


0.979910 -1.085444 -0.086820 


0.106525 -0.945238 




6 


0.970208 -1.078820 -0.086820 


0.082353 -0.957939 




7 


0.958651 -1.072610 -0.086820 


0.060678 -0.969186 




8 


0.945283 -1.066794 -0.086820 


0.040815 -0.979380 




9 


0.930159 -1.061323 -0.086820 


0.022307 -0.988783 




10 


0.913336 -1.056157 -0.086820 


0.004851 -0.997572 




11 


0.894883 -1.051250 -0.086820 


-0.011780 -1.005873 




12 


0.874870 -1.046545 -0.086820 


-0.027772 -1.013791 




13 


0.853379 -1.041994 -0.086820 


-0.043274 -1.021408 




14 


0.830493 -1.037547 -0.086820 


-0.058416 -1.028793 




15 


0.806302 -1.033154 -0.086820 


-0.073319 -1.036011 




16 


0.780903 -1.028773 -0.086820 


-0.088096 -1.043118 




17 


0.754395 -1.024362 -0.086820 


-0.102850 -1.050167 




18 


0.726883 -1.019887 -0.086820 


-0.117685 -1.057206 




19 


0.698476 -1.015315 -0.086820 


-0.132698 -1.064283 




20 


0.669285 -1.010628 -0.086820 


-0.147973 -1.071435 




21 


0.639427 -1.005803 -0.086820 


-0.163603 -1.078704 




22 


0.609018 -1.000836 -0.086820 


-0.179660 -1.086122 




23 


0.578179 -0.995714 -0.086820 


-0.196223 -1.093720 




24 


0.547031 -0.990441 -0.086820 


-0.213352 -1.101523 




25 


0.515698 -0.985017 -0.086820 


-0.231108 -1.109553 




26 


0.484302 -0.979449 -0.086820 


-0.249545 -1.117830 




27 


0.452969 -0.973745 -0.086820 


-0.268710 -1.126370 




28 


0.421821 -0.967910 -0.086820 


-0.288653 -1.135188 




29 


0.390982 -0.961956 -0.086820 


-0.309418 -1.144298 




30 


0.360573 -0.955883 -0.086820 


-0.331061 -1.153716 




31 


0.330715 -0.949692 -0.086820 


-0.353649 -1.163464 




32 


0.301524 -0.943373 -0.086820 


-0.377267 -1.173570 




33 


0.273117 -0.936907 -0.086820 


-0.402032 -1.184074 




34 


0.245605 -0.930250 -0.086820 


-0.428106 -1.195034 




35 


0.219097 -0.923337 -0.086820 


-0.455707 -1.206527 




36 


0.193698 -0.916066 -0.086820 


-0.485152 -1.218668 




37 


0.169507 -0.908273 -0.086820 


-0.516866 -1.231611 




38 


0.146621 -0.899722 -0.086820 


-0.551434 -1.245566 





Figure 2.3a UPOT output 
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39 


0.125130 


-C. 890054 -0.086820 


-0.589671 


-1.260822 






40 


0.105117 


-0.879738 -0.086820 


-0.632708 


-1.277774 






a: 


0.086664 


-0.864972 -0.086820 


-0 . 682121 


-1.296966 






A2 


0.069841 


-0.347530 -0.086820 


-0.740144 


-1.319145 






43 


0.054717 


-0.824472 -0.086820 


-0.810000 


-1.345363 






44 


0.041349 


-0.792611 -0.086820 


-0.896386 


-1.377093 






45 


0.029792 


-0.7<6398 -0.086820 


-1 .006151 


-1.416387 






46 


0.020090 


-0.675433 -0.086820 


-1 .148872 


-1.465903 






47 


0.012282 


-0.558323 -0.086820 


-1.334727 


-1.527981 






48 


0.006399 


-0.346398 -0.086820 


-1.554189 


-1.598183 






49 


0.002465 


0.0750C8 -0.086820 


-1.656311 


-1.629819 






50 


0.000493 


0.854855 -0.086820 


-0.999986 


-1.414209 






51 


0.000493 


1.561213 -0.086820 


0.449705 


-0.741819 






52 


0.002465 


1.669507 -0.086820 


0.997431 


-0.050689 






53 


0.006399 


1 .558290 -0.086820 


0.875814 


0.352401 






54 


0.012282 


1.445845 -0.086820 


0.674864 


0.570207 






55 


0.020090 


1.357269 -0.086820 


0.513069 


0.697804 






56 


0.029792 


1.287651 -0.086820 


0.392567 


0.779380 






57 


0.041349 


1.231267 -0.086820 


0.302348 


0.835256 






58 


0.054717 


1.184229 -0.086820 


0.233416 


0.875548 






59 


0.069841 


1 . 144055 -0.086820 


0.179673 


0.905719 






60 


0.086664 


1.109153 -0.086820 


0.137083 


0.928933 






61 


0.105117 


1 .078474 -0.086820 


0.102934 


0. 947136 






62 


0.12513C 


1.051316 -0.086820 


0.075364 


0.961580 






63 


0.146621 


1.027193 -0.086820 


0.053053 


0.973112 






64 


0.169507 


1.CC5756 -0.086820 


0.035030 


0.982329 






65 


0.193698 


0.986749 -0.086820 


0.020565 


0. 989664 






66 


0.219097 


0.969966 -0.086820 


0.009091 


0.995444 






67 


0.245605 


0.955241 -0.086820 


0.000149 


0.999926 






68 


0.273117 


0.942428 -0.086820 


-0.006640 


1.003314 






69 


0.301524 


0.931390 -0.086820 


-0.011593 


1.005780 






70 


0.330715 


0.921993 -0.086820 


-0.014984 


1.007464 






71 


0.360573 


0.914106 -0.086820 


-0.017049 


1.008488 






72 


0.390982 


0.907597 -0.086820 


-0.017994 


1.008957 






73 


0.421821 


0.902331 -0.086820 


-0.018001 


1.008960 






74 


0.452969 


0.898173 -0.086820 


-0.017221 


1.008574 






75 


0.484302 


0.894989 -0.086820 


-0.015784 


1.007861 






76 


0.515698 


0.892645 -0.086820 


-0.013805 


1 .006879 






7.7 


0.547031 


0.891016 -0.086820 


-0.011369 


1.005668 






78 


0.578179 


0.889982 -0.086820 


-0.008545 


1.004263 






79 


0.b09018 


0.889434 -0.086820 


-0.005387 


1.002690 






80 


0.639427 


0.889276 -0.086820 


-0.001926 


1.000962 






81 


0.669285 


0.889426 -0.086820 


0.001821 


0.999089 






S2 


0.698476 


0.889824 -0.086820 


0.005855 


0.997068 






83 


0.726883 


0.890423 -0.086820 


0.010190 


0. 994892 






84 


0.754395 


0.891196 -0.086820 


0.014851 


0.992547 






85 


0.780903 


0.892138 -0.086820 


0.019880 


0.990010 






86 


0.806302 


0.893264 -0.086820 


0.025328 


0.987255 






87 


0.830493 


0.894605 -0.086820 


0.031258 


0.984247 






88 


0.853379 


0.896208 -0.086820 


0.037745 


0.980946 






89 


0.874870 


0.898136 -0.086820 


0.044881 


0.977302 






90 


0.894883 


0.900467 -0.086820 


0.052775 


0.973255 






91 


0.913336 


0.903280 -0.086820 


0.061567 


0.968728 






92 


0.930159 


0.906675 -0.086820 


0.071417 


0.963630 






93 


0.945283 


0.910746 -0.086820 


0.082565 


0.957828 






94 


0.958651 


0.915603 -0.086820 


0.095332 


0.951141 






95 


0.970208 


0.921351 -0.086820 


0.110175 


0.943306 






96 


0. 979910 


0.928123 -0.086820 


0.127830 


0.933900 






97 


0.987718 


0.936086 -0.086820 


0.149558 


0.922194 






98 


0.993600 


0.945516 -0.086820 


0.177698 


0.906809 






99 


0.997535 


0.957020 -0.086820 


0.217400 


0.884647 






100 


0.999507 


0.971841 -0.086820 


0.280753 


0.848084 






... 


BEGIN UNSTEADY FLOW SOLUTION 


.... 








istep alpha 


time nitr 


cl , cd 


, cm 






1 


-3.0000 


0.0000 1 


-0.3479 


0.0002 


0.0403 




2 


-2.9836 


0.1540 0 


-0.3372 


-0.0005 


•0.0424 




3 


-2.9344 


0.3080 0 


-0.3247 


-0.0012 


0.0434 




4 


-2.8532 


0.4620 0 


-0.3101 


-0.0019 


0.0441 





Figure 2 . 3b UPOT output 
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stdin 



Page 3 



5 


- 2.7406 


0.6160 


c 


- 0.2935 


- 0.0027 


- 0.0444 


6 


- 2.5981 


0.7700 


c 


- 0.2751 


- 0.0033 


- 0.0444 


7 


- 2.4271 


0.9240 


0 


- 0.2549 


- 0.0039 


- 0 . C 441 


8 


- 2.2294 


1.0780 


r 


- 0.2332 


- 0.0044 


- 0.0434 


9 


- 2.0074 


1.2320 


c 


- 0.2101 


- 0.0047 


- 0.0424 


10 


-i . 7634 


1.3860 




- 0.1858 


- 0.0048 


- 0.0410 


11 


- 1 . 5000 


1.5400 




- 0.1605 


- 0.0048 


- 0.0393 


12 


- 1 .2202 


1.6940 




- 0.1345 


- 0.0046 


- 0.0373 


13 


- 0.9271 


1.8480 


0 


- 0.1080 


- 0.0043 


- 0.0350 


14 


- 0.6237 


2.0020 


c 


- 0.0812 


- 0.0038 


- 0.0323 


15 


- 0.3136 


2.1560 


0 


- 0.0545 


- 0.0031 


- 0.0295 


16 


0.0000 


2.3100 


; 


- 0.0280 


- 0.0023 


- 0.0264 


17 


0.3136 


2.4640 


c 


- 0.0021 


- 0.0014 


- 0.0231 


18 


0. 6237 


2.6180 


c 


0.0230 


- 0.0004 


- 0.0197 


19 


0.9271 


2.7720 


c 


0.0471 


0.0006 


- 0.0161 


20 


1.2202 


2.9260 


0 


0.0699 


0.0016 


- 0.0124 


21 


1.5000 


3.0800 




0.0911 


0.0025 


- 0.0087 


22 


1.7634 


3.2340 


c 


0.1107 


0.0034 


- 0.0050 


23 


2.0074 


3.3880 


n 


0. 1284 


0.0042 


- 0.0013 


24 


2 .2294 


3.5420 


0 


0,1440 


0.0048 


0.0024 


25 


2.4271 


3.6960 


c 


0.1573 


0.0053 


0.0059 


26 


2.5981 


3.8500 




0.1683 


0.0056 


0.0093 


27 


2.7406 


4.0040 


c 


0.1769 


0.0057 


0.0126 


28 


2.8532 


4.1580 


0 


0 . 1829 


0.0056 


0.0156 


29 


2.9344 


4.3120 


c 


0. 1863 


0.0054 


0.0184 


30 


2 . 9836 


4.4660 


c 


0. 1872 


0.0050 


0.0210 


31 


3.0000 


4.6200 


c 


0.1855 


0.0045 


0.0232 


32 


2.9836 


4.7740 


0 


0,1812 


0.0038 


0.0252 


33 


2.9344 


4.9280 


0 


0.1745 


0.0031 


0.0268 


34 


2.8532 


5.0820 


c 


0.1653 


0.0024 


0.0280 


35 


2.7406 


5.2360 


i 


0.1538 


0.0016 


0.0289 


36 


2.5981 


5.3900 


0 


0.1402 


0.0008 


0.0295 


37 


2.4270 


5.5440 


0 


0.1246 


0.0002 


0.0296 


38 


2 .2294 


5.6980 


0 


0.1071 


- 0.0004 


0.0294 


39 


2.0074 


5.8520 


c 


0.0880 


- 0.0009 


0.0288 


40 


1.7634 


6.0060 


c 


0.0676 


- 0.0012 


0.0278 


41 


1.5000 


6.1600 


G 


0.0459 


- 0.0014 


0.0265 


42 


1.2202 


6.3140 


0 


0.0234 


- 0.0014 


0.0249 


43 


0.9270 


6.4680 


0 


0.0001 


- 0.0013 


0.0229 


44 


0.6237 


6 . 62-20 


0 


- 0.0235 


- 0.0010 


0.0206 


45 


0.3136 


6.7760 


0 


- 0.0472 


- 0.0006 


0.0181 


46 


0.0000 


6.9300 




- 0.0709 


0.0000 


0.0153 


47 


- 0.3136 


7.0840 


0 


- 0.0941 


0.0006 


0.0123 


48 


- 0.6237 


7.2380 


0 


- 0.1166 


0.0014 


0.0092 


49 


- 0.9271 


7.3920 


0 


- 0.1382 


0.0021 


0.0059 


50 


- 1.2202 


7.5460 


0 


- 0.1586 


0.0028 


0.0025 


51 


- 1.5000 


7.7000 


1 


- 0.1776 


0.0035 


- 0.0010 


52 


- 1.7634 


7.8540 


0 


- 0.1950 


0.0041 


- 0.0045 


53 


- 2.0074 


8.0080 


0 


- 0.2105 


0.0047 


- 0.0080 


54 


- 2.2294 


8. 1620 


0 


- 0.2241 


0.0051 


- 0.0114 


55 


- 2.4271 


8.3160 


0 


- 0.2356 


0.0053 


- 0.0147 


56 


- 2.5981 


8.4700 


0 


- 0.2448 


0.0054 


- 0.0179 


57 


- 2.7406 


8.6240 


1 


- 0.2516 


0.0053 


- 0.0210 


58 


- 2. 8532 


8.7780 


0 


- 0.2559 


0.0051 


- 0.0239 


59 


- 2.9344 


8.9320 


0 


- 0.2578 


0.0047 


- 0.0265 


60 


- 2.9836 


9.0860 


0 


- 0.2571 


0,0042 


- 0.0289 


61 


- 3.0000 


9.2400 


0 


- 0.2539 


0.0035 


- 0.0309 


62 


- 2.9836 


9.3940 


0 


- 0.2482 


0.0028 


- 0.0327 


63 


- 2 . 9344 


9.5480 


0 


- 0.2400 


0.0020 


- 0.0342 


64 


- 2.8532 


9.7020 


0 


- 0.2295 


0.0011 


- 0.0353 


65 


- 2.7406 


9.8560 


0 


- 0.2168 


0.0003 


- 0.0360 


66 


- 2.5981 


10.0100 


0 


- 0.2019 


- 0.0005 


- 0.0364 


67 


- 2.4271 


10.1640 


0 


- 0.1851 


- 0.0012 


- 0.0365 


68 


- 2.2294 


10.3180 


0 


- 0.1664 


- 0.0017 


- 0.0361 


69 


- 2.0074 


10.4720 


0 


- 0.1463 


- 0.0022 


- 0.0354 


70 


- 1.7634 


10.6260 


0 


- 0.1247 


- 0.0025 


- 0.0343 


71 


- 1 . 5000 


10.7800 


0 


- 0.1020 


- 0.0026 


- 0.0329 


72 


- i .2202 


10.9340 


0 


- 0.0784 


- 0.0026 


- 0.0311 


73 


- 0.9271 


11.0880 


0 


- 0.0542 


- 0.0024 


- 0.0290 


74 


- 0.6237 


11.2420 


0 


- 0.0296 


- 0.0020 


- 0.0267 


75 


- 0.3136 


11.3960 


0 


- 0.0050 


- 0.0015 


- 0.0240 


76 


0.0000 


11.5500 


0 


0.0195 


- 0.0008 


- 0.0211 


77 


0.3136 


11.7040 


0 


0.0436 


- 0.0001 


- 0.0181 



Figure 2.3c UPOT output 
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stdin 


Page 4 


78 


0.6237 


11.8580 


0 


0.0670 0.0007 -0.0148 




^9 


C.9270 


12,0120 


0 


0.0894 0.0016 -0.0114 




8C 


1 .2202 


12.1660 


1 


0.1106 0.0024 -0.0079 




81 


: .5303 


12.3200 


0 


0.1304 0.0033 -0.0044 




32 


:.7633 


12.4740 


0 


0.1486 0.0040 -0.0008 




83 


2. 3074 


12,6280 


0 


0.1649 0,0046 0.0028 




84 


2.2294 


12 .7820 


0 


0.1792 0.0051 0.0063 




85 


2.4270 


12.9360 


1 


0.1913 0.0055 0.0097 




86 


2.5981 


13.0900 


0 


0.2012 0.0057 0.0130 




87 


2.7406 


13.2440 


0 


0.2087 0.0057 0.0161 




88 


2.8532 


13 . 3980 


0 


0.2137 0.0056 0.0191 




89 


2.9344 


13.5520 


0 


0.2161 0.0053 0.0218 




9C 


2.9836 


13.7060 


0 


0.2160 0.0048 0.0242 




91 


3.0300 


13.8600 


0 


0.2134 0.0042 0.0263 




92 


2.9336 


14.0140 


1 


0.2083 0.0036 0.0282 




93 


2.9344 


14.1680 


0 


0.2007 0.0028 0.0297 




94 


2.8532 


14.3220 


0 


0.1907 0.0020 0.0309 




9: 


2. M06 


14.4760 


0 


0.1785 0.0012 0.0317 




96 


2.5981 


14 . 6300 


0 


0.1641 0.0005 0.0321 




97 


2.4271 


14.7840 


0 


0.1478 -0.0002 0.0322 




98 


2.2294 


14.9380 


0 


0.1297 -0.0008 0.0319 




99 


2.0074 


15.0920 


0 


0.1100 -0.0013 0.0312 




IOC 


1.7634 


15.2460 


0 


0.0889 -0.0016 0.0302 




ICl 


1.5000 


15.4000 


0 


0.0666 -0.0017 0.0288 




102 


1 .2202 


15.5540 


0 


0.0435 -0.0017 0.0271 




103 


0.9271 


15,7080 


0 


0.0197 -0.0016 0.0251 




104 


3.6238 


15.8620 


0 


-0.0044 -0.0013 0.0228 




105 


0.3136 


16.0160 


0 


-0.0287 -0.0008 0.0202 




106 


3.0000 


16.1700 


0 


-0.0528 -0.0002 0.0173 




107 


-0.3136 


16.3239 


1 


-0.0765 0.0005 0.0143 




109 


-0.6237 


16.4779 


0 


-0.0994 0.0012 0.0111 




109 


-0.9270 


16.6319 


0 


-0.1215 0.0020 0.0078 




no 


-1.2202 


16.7859 


0 


-0.1423 0.0028 0.0043 




111 


-1.5000 


16.9399 


0 


-0.1617 0.0035 0.0008 




112 


-1.7633 


17,0939 


1 


-0.1795 0.0042 -0.0027 




113 


-2.0074 


17.2479 


0 


-0.1955 0.0047 -0.0063 




114 


-2 .2294 


17.4019 


0 


-0.2095 0.0052 -0.0097 




115 


-2.4270 


17.5559 


0 


-0.2212 0.0054 -0.0131 




116 


-2.5981 


17.7099 


0 


-0.2308 0.0056 -0.0164 




117 


-2.7406 


17.8639 


0 


-0.2379 0.0055 -0.0195 




118 


-2.8532 


18.0179 


1 


-0.2426 0.0053 -0.0224 




119 


-2.9344 


18.1719 


0 


-0.2447 0,0049 -0.0250 




120 


-2.9836 


18,3259 


0 


-0,2444 0.0044 -0.0274 








PHASE SHIFT 


ANALYSIS 








FREO = 


0. 


6800000 




AMPLI 


T'JDE; clamp, cmamp 




0.2304234 3 . 4 331881 E-02 




PHASE 


; clp, cmp 




184.9092 -37.54200 




AVERAGE DRAG, TOTAL DRAG 




1 .5421067E-03 9 . 406851 2E-02 




ETAS, 


WBAR 






-0.2168084 -7 . 1127615E-03 




EREO SWEEP 










FREQ = 


0.690000 








STEADY FLOW SOLUTION AT ALPHA = -3.000000 




1 


0.999507 


-1.116717 


-0. 


086820 0.280753 -0.848084 




2 


0.997535 


-1 . 108171 


-0. 


086820 0.213074 -0.887088 




3 


0.993600 


-1.100064 


-0. 


086820 0.168200 -0,912031 




4 


0.987718 


-1.092523 


-0, 


086820 0.134381 -0.930386 




5 


0.979910 


-1.085444 


-0. 


086820 0.106525 -0.945238 




6 


0.970208 


-1.078820 


-0. 


086820 0.082353 -0.957939 




7 


0.958651 


-1.072610 


-0. 


086820 0.060678 -0.969186 




8 


0.945283 


-1.066794 


-0. 


086820 0.040815 -0.979380 




9 


0.930159 


-1.061323 


-0. 


086820 0.022307 -0.988783 




10 


0.913336 


-1.056157 


-0. 


086820 0.004851 -0.997572 




11 


0.894883 


-1.051250 


-0. 


086820 -0.011780 -1.005873 




12 


0.874870 


-1 .046545 


-0. 


086820 -0.027772 -1.013791 





Figure 2.3d UPOT output 
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stdin 


Page 5 


13 


0.853379 


- 1.041994 - 0.086820 


- 0.043274 


- 1.021408 




14 


0.830493 


- 1 . C 37547 - 0.086820 


- 0.058416 


- 1.028793 




15 


0.806302 


- 1.033154 - 0.086820 


- 0.073319 


- 1 .036011 




16 


0.780903 


- 1 . C 28773 - 0.086820 


- 0.088096 


- 1.043118 




17 


0.754395 


- 1.024362 - 0.086820 


- 0.102850 


- 1.050167 




18 


0.726883 


- 1.019887 - 0.086820 


- 0.117685 


- 1.057206 




19 


0.698476 


- 1.015315 - 0.086820 


- 0.132698 


- 1 .064283 




20 


0.669285 


- 1.010628 - 0.086820 


- 0.147973 


- 1.071435 




21 


0.639427 


- 1.005803 - 0.086820 


- 0.163603 


- 1.078704 




22 


0.609018 


- 1.000836 - 0.086820 


- 0.179660 


- 1,086122 




23 


0.578179 


- 0.995714 - 0.086820 


- 0.196223 


- 1.093720 




24 


0.547031 


- 0.990441 - 0.086820 


- 0.213352 


- 1 ,101523 




25 


0.515698 


- 0.985017 - 0.086820 


- 0.231108 


- 1.109553 




26 


0.484302 


- 0.979449 - 0,086820 


- 0.249545 


- 1.117830 




27 


0.452969 


- 0.973745 - C . 086820 


- 0.268710 


- 1 .126370 




28 


0.421821 


- 0.967910 - 0.086820 


- 0.288653 


- 1 .135188 




29 


0. 390982 


- 0.961956 - 0.086820 


- 0.309418 


- 1.144298 




30 


0.360573 


- 0.955883 - C . 086820 


- 0.331061 


- 1.153716 




31 


0.330715 


- 0.949692 - 0.086820 


- 0.353649 


- 1.163464 




32 


0.301524 


- 0.943373 - 0.086820 


- 0.377267 


- 1 .1 73570 




33 


0.273117 


- 0.936907 - 0.086820 


- 0.402032 


- 1.184074 




34 


0.245605 


- 0.930250 - 0.086820 


- C . 428106 


- 1.1 95034 




35 


0.219097 


- 0.923337 - 0.086820 


- 0.455707 


- 1.206527 




36 


0.193698 


- 0.916066 - 0.086820 


- 0.485152 


- 1,218668 




37 


0.169507 


- 0.908273 - 0.086820 


- 0.516866 


- 1.231611 




38 


0.146621 


- 0.899722 - 0.086820 


- 0.551434 


- 1.245566 




39 


0.125130 


- 0.890054 - 0.086820 


- 0,589671 


- 1.260822 




40 


0.105117 


- 0.878738 - 0.086820 


- 0.632708 


- 1.277774 




41 


0.086664 


- 0.864972 - 0.086820 


- 0.682121 


- 1.296966 




42 


0.069841 


- 0.847530 - 0.086820 


- 0.740144 


- 1.319145 




43 


0.054717 


- 0.824472 - 0.086820 


- 0.810000 


- 1.345363 




44 


0.041349 


- 0.792611 - 0.086820 


- 0.896386 


- 1.377093 




45 


0.029792 


- 0.746398 - 0.086820 


- 1.006151 


- 1.416387 




46 


0.020090 


- 0.675433 - 0.086820 


- 1.148872 


- 1.465903 




47 


0.012282 


- 0.558323 - 0.086820 


- 1.334727 


- 1.527981 




48 


0.006399 


- 0.346398 - 0.086820 


- 1.554189 


- 1.598183 




49 


0.002465 


0.075008 - 0.086820 


- 1.656311 


- 1.629819 




50 


0.000493 


0.854855 - 0.086820 


- 0.999986 


- 1.414209 




51 


0.000493 


1.561213 - 0.086820 


0.449705 


- 0,741819 




52 


0.002465 


1.669507 - 0 . 08^820 


0.997431 


- 0.050689 




53 


0.006399 


1.558290 - 0.086820 


0.875814 


0.352401 




54 


0.012282 


1.445845 - 0.086820 


0.674864 


0.570207 




55 


0.020090 


1.357269 - 0.086820 


0.513069 


0.697804 




56 


0.029792 


1.287651 - 0.086820 


0.392567 


0,779380 




57 


0.041349 


1.231267 - 0.086820 


0.302348 


0.835256 




58 


0.054717 


1.184229 - 0.086820 


0.233416 


0.875548 




59 


0.069841 


1.144055 - 0.086820 


0.179673 


0.905719 




60 


0.086664 


1.109153 - 0.086820 


0.137083 


0.928933 




61 


0.105117 


1.078474 - 0.086820 


0.102934 


0.947136 




62 


0.125130 


1.051316 - 0.086820 


0.075364 


0.961580 




63 


0.146621 


1.027193 - 0.086820 


0.053053 


0.973112 




64 


0.169507 


1.005756 - 0.086820 


0.035030 


0.982329 




65 


0.193698 


0.986749 - 0.086820 


0.020565 


0.989664 




66 


0.219097 


0.969966 - 0.086820 


0.009091 


0.995444 




67 


0.245605 


0.955241 - 0.086820 


0.000149 


0.999926 




68 


0.273117 


0.942428 - 0.086820 


- 0.006640 


1.003314 




69 


0.301524 


0.931390 - 0.086820 


- 0.011593 


1.005780 




70 


0.330715 


0.921993 - 0.086820 


- 0.014984 


1.007464 




71 


0.360573 


0.914106 - 0.086820 


- 0.017049 


1.008488 




72 


0.390982 


0.907597 - 0.086820 


- 0.017994 


1.008957 




73 


0.421821 


0.902331 - 0.086820 


- 0.018001 


1.008960 




74 


0.452969 


0.898173 - 0.086820 


- 0.017221 


1.008574 




75 


0.484302 


0.894989 - 0.086820 


- 0.015784 


1.007861 




76 


0.515698 


0.892645 - 0.086820 


- 0.013805 


1.006879 




77 


0.547031 


0.891016 - 0.086820 


- 0.011369 


1.005668 




78 


0.578179 


0.889982 - 0.086820 


- 0.008545 


1.004263 




79 


0.609018 


0.889434 - 0.086820 


- 0.005387 


1.002690 




80 


0.639427 


0.889276 - 0.086820 


- 0.001926 


1.000962 




81 


0.669285 


0.889426 - 0.086820 


0.001821 


0.999089 




82 


0.698476 


0.889824 - 0.086820 


0.005855 


0.997068 




83 


0.726883 


0.890423 - 0.086820 


0.010190 


0. 994892 




84 


0.754395 


0.891196 - 0.086820 


0.014851 


0.992547 




85 


0.780903 


0.892138 - 0.086820 


0.019880 


0.990010 





Figure 2.3e UPOT output 
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stdin 
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86 


0.806302 


0.893264 


- 0.086820 


0.025328 


0.987255 






87 


0.830493 


0.894605 


- 0.086820 


0.031258 


0.984247 






88 


0.853379 


0.896208 


- 0.086820 


0.037745 


0.980946 






89 


C . 874870 


0.898136 


- 0.086820 


0.044881 


0.977302 






90 


0.894883 


0.900467 


- 0.086820 


0.052775 


0.973255 






9 : 


0.913336 


0.903280 


- 0.086820 


0.061567 


0.968728 






92 


0.930159 


0.906675 


- 0.086820 


0.071417 


0.963630 






93 


0.945283 


0.910746 


- 0.086820 


0.082565 


0.957828 






94 


0.958651 


0.915603 


- 0.086820 


0.095332 


0.951141 






95 


0.970208 


0.921351 


- 0.086820 


0.110175 


0.943306 






96 


0.979910 


0.928123 


- 0.086820 


0.127830 


0.933900 






97 


0.987718 


0.936086 


- 0.086820 


0.149558 


0.922194 






98 


0.993600 


0.945516 


- 0.086820 


0.177698 


0.906809 






99 


0.997535 


0.957020 


- 0.086820 


0.217400 


0.884647 






100 


0.999507 


0.971841 


- 0.086820 


0.280753 


0.848084 






... 


BEGIN UNSTEADY PLOW 


SOLUTION 


.... 








isteo alpha 


Line 


nit r 


cl , cd 


, cm 








- 3.0000 


0.0000 


1 


- 0.3479 


0.0002 


- 0.0403 




2 


- 2.9836 


0.1518 


0 


- 0.3370 


- 0.0005 


- 0.0425 




3 


- 2.9344 


0.3035 


0 


- 0.3243 


- 0.0012 


- 0.0435 




4 


- 2.8532 


0.4553 


0 


- 0.3096 


- 0.0020 


- 0.0442 




5 


- 2.7406 


0,6071 


0 


- 0.2929 


- 0.0027 


- 0.0445 




6 


- 2.5981 


0.7588 


0 


- 0.2744 


- 0.0034 


- 0.0446 




7 


- 2.4271 


0.9106 


0 


- 0.2542 


- 0.0039 


- 0.0442 




8 


- 2.2294 


1.0624 


0 


- 0.2323 


- 0.0044 


- 0.0436 




9 


- 2.0074 


1.2141 


0 


- 0.2091 


- 0.0047 


- 0.0426 




10 


- 1.7634 


1.3659 


0 


- 0.1848 


- 0.0049 


- 0.0412 




i : 


- 1.5000 


1.5177 


0 


- 0.1595 


- 0.0049 


- 0.0395 




12 


- 1.2202 


1.6694 


0 


- 0.1334 


- 0.0047 


- 0.0375 




13 


- 0.9271 


1.8212 


0 


- 0.1069 


- 0.0043 


- 0.0352 




14 


- 0.6237 


1.9730 


0 


- 0.0802 


- 0.0038 


- 0.0326 




15 


- 0.3136 


2.1247 


0 


- 0.0535 


- 0.0031 


- 0.0297 




16 


0.0000 


2.2765 


1 


- 0.0271 


- 0.0023 


- 0.0267 




17 


0.3136 


2.4283 


0 


- 0.0012 


- 0.0014 


- 0.0234 




18 


0.6237 


2.5801 


0 


0.0238 


- 0.0004 


- 0.0199 




19 


0.9271 


2.7318 


0 


0.0478 


0.0006 


- 0.0164 




20 


1.2202 


2.8836 


0 


0.0704 


0.0016 


- 0.0127 




21 


1.5000 


3.0354 


1 


0.0916 


0.0026 


- 0.0090 




22 


1.7634 


3.1871 


0 


0.1110 


0.0034 


- 0.0052 




23 


2.0074 


3.3389 


0 


0.1285 


0.0042 


- 0.0015 




24 


2.2294 


3.4907 


0 


0.1439 


0.0048 


0.0022 




25 


2.4271 


3.6424 


0 


0.1571 


0.0053 


0.0057 




26 


2.5981 


3.7942 


1 


0.1679 


0.0056 


0.0092 




27 


2.7406 


3. 9460 


0 


0.1763 


0.0057 


0.0124 




28 


2.8532 


4.0977 


0 


0.1822 


0.0056 


0.0155 




29 


2.9344 


4.2495 


0 


0.1855 


0.0054 


0.0183 




30 


2.9836 


4.4013 


0 


0.1862 


0.0050 


0.0209 




31 


3.0000 


4.5530 


0 


0.1843 


0.0045 


0.0231 




32 


2.9836 


4.7048 


0 


0.1799 


0.0038 


0.0251 




33 


2.9344 


4.8566 


0 


0.1729 


0.0031 


0.0267 




34 


2.8532 


5.0083 


0 


0.1636 


0.0023 


0.0280 




35 


2.7406 


5.1601 


1 


0.1520 


0.0016 


0.0289 




36 


2.5981 


5.3119 


0 


0.1383 


0.0008 


0.0295 




37 


2.4271 


5.4636 


0 


0.1226 


0.0001 


0.0296 




38 


2.2294 


5.6154 


0 


0.1050 


- 0.0005 


0.0294 




39 


2.0074 


5.7672 


0 


0.0859 


- 0.0009 


0.0288 




40 


1.7634 


5.9189 


0 


0.0654 


- 0.0013 


0.0279 




41 


1.5000 


6.0707 


0 


0.0437 


- 0.0014 


0.0266 




42 


1 .2202 


6.2225 


0 


0.0211 


- 0.0014 


0.0250 




43 


0.9271 


6.3742 


0 


- 0.0021 


- 0.0013 


0.0230 




44 


0.6237 


6.5260 


0 


- 0.0257 


- 0.0010 


0.0207 




45 


0.3136 


6.6778 


0 


- 0.0494 


- 0.0006 


0.0182 




46 


0.0000 


6.8295 


1 


- 0.0729 


0.0000 


0.0154 




47 


- 0.3136 


6.9813 


0 


- 0.0960 


0.0007 


0.0125 




48 


- 0.6237 


7.1331 


0 


- 0.1185 


0.0014 


0.0093 




49 


- 0.9270 


7.2849 


0 


- 0.1399 


0.0021 


0.0060 




50 


-1 .2202 


7.4366 


0 


- 0.1602 


0.0029 


0.0026 




51 


- 1.5000 


7.5884 


1 


- 0.1791 


0.0036 


- 0.0009 





Figure 2.3f UPOT output 
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stdin 



52 


< 1.7634 


7.7402 


0 


- 0.1963 


0.0042 


- 0.0044 


53 


- 2.0074 


7.8919 


0 


- 0.2117 


0.0047 


- 0.0079 


54 


- 2.2294 


8.0437 


0 


- 0.2251 


0.0051 


- 0.0113 


55 


- 2 . 427 C 


8.1955 


0 


- 0.2364 


0.0053 


- 0.0147 


56 


- 2.5981 


8.3472 


1 


- 0.2454 


0.0054 


- 0.0179 


57 


- 2.7406 


8.4990 


C 


- 0.2520 


0 . 0 C 53 


- 0.0210 


58 


- 2.8532 


8.6508 


0 


- 0.2562 


0.0051 


- 0.0238 


59 


- 2.9344 


8.8025 


c 


- 0.2578 


0.0047 


- 0.0265 


60 


- 2.9836 


8.9543 


0 


- 0.2570 


0.0042 


- 0.0289 


61 


- 3.0000 


9.1061 


0 


- 0.2536 


0.0035 


- 0.0310 


62 


- 2.9836 


9.2578 


0 


- 0.2477 


0.0028 


- 0.0328 


63 


- 2.9344 


9.4096 


0 


- 0.2394 


0.0019 


- 0.0342 


64 


- 2.8532 


9.5614 


0 


- 0.2288 


0.0011 


- 0.0354 


65 


- 2.7406 


9.7131 


c 


- 0.2159 


0.0003 


- 0.0361 


66 


- 2.5981 


9.8649 


0 


- 0.2009 


- C .0005 


- 0.0365 


67 


- 2 . 427 : 


10.0167 


0 


- 0.1839 


- 0.0012 


- 0.0366 


68 


- 2.2294 


10.1684 


0 


- 0.1652 


- 0.0018 


- 0.0362 


69 


- 2.0074 


10. 3202 


V 


- 0.1450 


- 0.0022 


- 0.0355 


70 


- 1.7634 


10.4720 




- 0.1233 


- 0.0025 


- 0.0345 


71 


- 1.5000 


10.6237 


c 


- 0.1006 


- 0.0027 


- C . 033 C 


72 


- 1.2202 


10.7755 


0 


- 0.0770 


- 0.0026 


- 0.0313 


73 


- 0.9271 


: 0. 9273 


0 


- 0.0528 


- 0.0024 


- 0.0292 


74 


- 0.6237 


11.0790 


0 


- 0.0282 


- 0.0020 


- 0.0269 


75 


- 0.3136 


11 .2308 


0 


- 0.0036 


- 0.0015 


- 0.0242 


76 


O.OOOC 


1 1 . 3826 


c 


0.0209 


- 0.0008 


- 0.0215 


77 


0.3136 


11.5343 


c 


0.0448 


- 0.0001 


- 0.0183 


78 


0.6237 


11.6861 


0 


0.0681 


0.0007 


- 0.0150 


79 


0.9270 


11.8379 


0 


0.0904 


0.0016 


- 0.0116 


80 


1.2202 


11.9897 


0 


0.1115 


0.0025 


- 0.0081 


81 


1.5000 


12.1414 


1 


0.1311 


0.0033 


- 0.0045 


82 


1.7634 


12.2932 


0 


0.1493 


0.0040 


- 0.0010 


83 


2.0074 


12.4450 


0 


0.1654 


0.0047 


0.0026 


84 


2.2294 


12.5967 


0 


0.1796 


0.0052 


0.0061 


85 


2.4270 


12.7485 


0 


0.1915 


0.0055 


0.0095 


86 


2.5981 


12.9003 


1 


0.2012 


0.0057 


0.0129 


87 


2.7406 


13.0520 


0 


0.2085 


0.0058 


0.0160 


88 


2.8532 


13.2038 


0 


0.2133 


0.0056 


0.0189 


89 


2.9344 


13.3556 


0 


0.2156 


0.0053 


0.0217 


90 


2.9836 


13.5073 


0 


0.2153 


0.0048 


0.0241 


91 


3.0000 


13.6591 


0 


0.2125 


0.0042 


0.0263 


92 


2.9836 


13.8109 


0 


0.2072 


0.0036 


0.0282 


93 


2.9344 


13.9626 


0 


0.1995 


0.0028 


0.0297 


94 


2.8532 


14.1144 


0 


0.1894 


0.0020 


0.0309 


95 


2.7406 


14.2662 


0 


0.1770 


0.0012 


0.0317 


96 


2.5981 


14.4179 


1 


0.1625 


0.0005 


0.0322 


97 


2.4271 


14.5697 


0 


0.1461 


- 0.0002 


0.0323 


98 


2.2294 


14.7215 


0 


0.1279 


- 0.0008 


0.0320 


99 


2.0074 


14 . 8.732 


0 


0.1081 


- 0.0013 


0.0313 


100 


1.7634 


15.0250 


0 


0.0870 


- 0.0016 


0.0303 


101 


1.5000 


15.1768 


0 


0.0647 


- 0.0018 


0.0289 


102 


1.2202 


15. 3285 


0 


0.0415 


- 0.0018 


0.0272 


103 


0.9271 


15.4803 


0 


0.0178 


- 0.0016 


0.0252 


104 


0.6237 


15.6321 


0 


- 0.0063 


- 0.0013 


0.0229 


105 


0.3136 


15.7838 


0 


- 0.0305 


- 0.0008 


0.0203 


106 


0.0000 


15.9356 


1 


- 0.0546 


- 0.0002 


0.0175 


107 


- 0.3136 


16.0874 


0 


- 0.0782 


0.0005 


0.0145 


108 


- 0.6237 


16.2391 


0 


- 0.1011 


0.0012 


0.0113 


109 


- 0.9270 


16.3909 


0 


- 0.1230 


0.0020 


0.0079 


no 


- 1.2202 


16.5427 


0 


- 0.1437 


0.0028 


0.0044 


111 


- 1.5000 


16.6945 


1 


- 0.1630 


0.0035 


0.0009 


112 


- 1.7634 


16.8462 


0 


- 0.1807 


0.0042 


- 0.0026 


113 


- 2.0074 


16.9980 


0 


- 0.1965 


0.0048 


- 0.0061 


114 


- 2.2294 


17.1498 


0 


- 0.2103 


0.0052 


- 0.0096 


115 


- 2.4270 


17.3015 


0 


- 0.2219 


0.0055 


- 0.0130 


116 


- 2.5981 


17.4533 


1 


- 0.2312 


0.0056 


- 0.0163 


117 


- 2.7406 


17 . 6051 


0 


- 0.2382 


0.0055 


- 0.0194 


118 


- 2.8532 


17.7568 


0 


- 0.2427 


0.0053 


- 0.0223 


119 


- 2.9344 


17.9086 


0 


- 0.2446 


0.0049 


- 0.0250 


120 


- 2.9836 


18.0604 


0 


- 0.2441 


0.0044 


- 0.0274 



PHASE SHIFT ANALYSIS 
FREQ = 0.6900000 



Figure 2.3g UPOT output 
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2. Difference Between UPOT and U2DIIF/Phase 



The input file format was changed along with the 
following : 



• The program can now analyze a pitch, plunge or ramp motion 
that starts from any minimum value of Alpha or plunge 
(h/2b) . Previously, the program only accepted the initial 
position of zero. This program does not need to go 
through the origin. 

• The phase portion of the program was changed to curve fit 
Cl and Cj^ to a cosine function: 

F( t) =Amp*Cos ((i) t+4>) (2.2) 

where. Amp = amplitude of motion, w = frequency of motion, 
<t> = phase angle between motion and the aerodynamic 
forces. This was done since the alpha and plunge values 
were allowed to start from a new zero position. 

• The phase portion uses the middle 180 degrees of the final 
360 degree cycle specified in the UPOT. in file. This 
change was done to capture an all positive area of the 
cosine curve for phasing analysis. The program integrates 
this portion of the cosine curve, and for proper code 
opera.tion the area under the curve must be kept to one 
sign. If the areas of integration were chosen to include 
both sides of the axis, then the code would produce errors 
near 90 and 270 degrees. 

3. UPOT Verification 

The code UPOT did not incorporate any drastic changes 
to the prior codes, but the original code had never been 
extensively compared to prior theories over a wide range of 
reduced frequencies. When conducting these comparisons, it is 
easy to become confused. This section will go through the 
comparisons slowly to help alleviate that problem. 
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^jjanel ^Theodorson^^t^ • 

The equation for reduced frequency is: 



, _ <j>2b „ _ (j>b 

— — — — ^ 2 \ — — ■ ■- 

P u ^ u 



(2.3) 



where: w = frequency of oscillation (rad/sec) 

b = half chord (units to match U) 

U = free stream velocity (units to match b) . 

The difference between Kp and lies in the fact that Kp 

calls for the full chord and calls for the half chord, 

hence, it is important to remember that Kp is twice K^. 
b. Aerodynamic Forces 

The aerodynamic forces problem of simple harmonic 
motion about an equilibrium position was solved theoretically 
by Theodorsen in NACA TR-496 [ref. 10] and outlined by Fung in 
[ref. 5]. The complex equations were simplified using the 
simple harmonic motion equation and resulted in the following: 

L=npb^u>^ ( A+a) LJ a+ [Lp- (c-e) L,] p) (2.4) 

M=7ipb"o)2( [w^-(l+a)L^] ^ + 

[M^-(l^a) (L„+M^) +(^+a)2Lj a (2.5) 

+ [Wp-(l+a)Lp-(c-e)M^+(c-e) i^^a)L^] p) 

L, M are the lift and moment per unit span of the 
airfoil about the elastic axis, b, h/b, a and a (radians) , are 



18 



shown in Figure (2.1)., L^, L^, Lp, and are defined by 

Scanlan [ref . 6 , pp . 412 - 424] for various values of and e. 
This analysis will not cover airfoil aileron combinations. 
Therefore ^ becomes zero and equations 2.4 and 2.5 reduce to: 

L=npb^(o^ ( 2 . 6 ) 



M=Jtpb«a>2( [M^-(l+a)Lj ^ + 
[M - ( ^ +a) (L +MJ + ( 1 +a) ^Lj a) e" 



(2.7) 



The UPOT panel code used the following equations 
in defining lift and moment: 

«|,^=tan-iii (2.8) 

^ =slRl*lle^ <j)^=tan-^4^ (2.9) 

4gb2 

where Rl and are the real and imaginary parts of Cj^, and Rj^ 
and Ijj are the real and imaginary parts of Cjj. For the same 
conditions the lift (L) and moment (M) should be the same for 
both the panel code and Theodorsen. This fact allows for 
comparison of the magnitude, real, imaginary and phase of lift 
and moment . 

For lift: (eqn 2.6) equals Lp (eqn 2.8) 
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After canceling e 



i (ei)t+0D . 



Ttpjb^o)^ [Lf^{h/b) + (L„- (Vi+a) L^) a] =2gi>jRl+ll (2.11) 

C^=^Rl^ll 



For pitch case, h/b = 0, 2.11 reduces to: 

7ipb^o)2 f (1/2+3) ^^) „] =2gbC„ (2.12) 

Substitute K^. = bo)/U for w* and q = into 

equation 2.12 which gives: 

2 Tc gbKt ( L„ - ( V2 + a ) L^) a =2 gbq„ (2.13) 

After cancelling and substituting Kp for K^.: 

^ [L„-(V 2 +a)Lja=Q„ (2-14) 



This relationship can be further broken down into the real and 
imaginary parts : 



Imag: [iL„- (Va+a) iL^] a =Q„sin ((J)^) (2.15) 



Plunge Case: O!=0 using equation 2.11 gives: 
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Real : 



( 2 . 16 ) 



nK, 



^ [L„ -(V 2 +a) LJ a=C._cos (<t>.) 



npb^(O^L^-^=2gbC^ 



( 2 . 17 ) 



The panel code uses h/2b for analysis because it uses full 
chord vice half chord. Therefore equation 2.17 becomes; 



27ipb^o)2L^(_) =2gbC^ 



( 2 . 18 ) 



Substituting as before for oj: 



2ngbKlL^2(-^) =2gbC^ 



Cancel and substitute Kp! 



(A)l =c 
^ 2 ^ ^ 2b^ ^ “ 



( 2 . 19 ) 



( 2 . 20 ) 



This can also be broken up into imaginary and real parts as 
before. 

MOMENT: 

Equating equations 2.7 and 2.9 results in: 



npb"a)2 ( (V 2 +a) LJ [M„- (V 2 +a) + (V 2 +a) a) 






For pitch: h/b = 0 



21 



( 2 . 22 ) 






^Ma '*'^1 



resulting in: 



C 



( [M„- (V 2 +a) (L^+Mfj) + (Va+a) a) =4gb^C^ 



Ma 



(2 



After substituting and cancelling: 



anK, 



^ [M„- (V 2 +a) (L„+M^) + (V 2 +a) 2 lJ =C^„ <2 



REAL Mv, = M 



anK, 



8 



^ [M„- (V 2 +a) (L„-^V 2 ) + (V 2 +a) ^Lj =Q„cos (<|)„„) <2 



IMAG: 



Mh = 0 



anK, 



[iM„-(V 2 +a) (iL„) + (V 2 +a) 2iL^] =Cj^„sin (4>^„) <2 



For plunge, a = 0, equation 21 reduces to : 
REAL: 






(2 



. 23 ) 



. 24 ) 



. 25 ) 



26 ) 



. 27 ) 



. 28 ) 
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0 



IMAG: Mh = 

(%ta)Lj=C^sin(4>^) (2.29) 

Comparisons are shown for various cases of pitch and plunge. 
The tables include pitch values of 1 (Tables 2. 2 -2. 3) and 6.7 
degrees (Tables 2. 4-2. 5), plunge (h/2b) values of .01 (Table 
2.8) and .0833 (Table 2. 6-2. 7). The graphs include 1 degree 
pitch (Figures 2.7-12, 2.13-20), 6.7 degree pitch (Figures 
2.21-28,2.29-36), . 01 h/2b plunge (Figures 2 . 53-56) , and .0833 
h/2b plunge (Figures 2.37-44, 2.45-52). 

4. Results 

The tables and graphs show that the panel code 
predicts the Theodorsen results accurately. An initial 
question that was first addressed for the comparisons was how 
many cycles to use for good consistent phase results. Initial 
runs were made at several different cycle values and the 
results are shown in Table 2.1. It can be seen that the panel 
code predicts Theodorsen' s results, using a cycle number of 
three. Increasing the cycle number just talces more computer 
time and only results in marginal increases in accuracy. 

The most glaring difference appears in the Ij^j c^jj^ 
comparisons of Figures 2.47-48. The panel code drops off 
sharply at the higher end of the reduced frequency spectrum. 
The reason for this is believed to be due to the magnitude of 
h/2b chosen for the comparison. The code was rerun for a 
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Comparison of Phase Calculations Using Various Cycles. 






(Pitch, 6.7 deg., NACA 0007, .37c, 50 panels top and bottom) 


















Kp 


# Cycles 


Cl Phase Angle 


Cm Phase Angle 


Cl Amplitude 


Cm amplituc 


e 


1.00 


1 


182.0537 


■54.409 


0.4884 


0.08937 




1.00 


2 


208.1592 


46.424 


0.5109224 


0.083932 




1.00 


3 


206.001 


-44.313 


0.51527 


0.083169 




% Dili. 213 




1.05% 


4.76% 


0.84% 


0.92% 


















1.00 


4 


204.9365 


-43.33887 


0.51668 


0.082951 




1.00 


5 


204.3955 


-42.86817 


0.5172907 


0.082865 1 


i 


1.00 


6 


204.0596 


-42.58497 


0.517598 


0.082822 




1.00 


7 


203.8291 


■42.3916 


0.517776 


0.082792 




1.00 


8 


203.6031 


■42.25684 


0.5 i 789 


0.082784 




% Dili. 7/8 




o.ii% 


0.32% 


0.02% 


0.01% 


















3.60 


1 2 


264.44052 


•59.79002 


1.0931 


0.212011 




3.60 


i 3 


261.9737 


•58.29977 


1.09954 


0.21936 




% OKI. 2/3 




0.94% 


2.56% 


0.59% 


3.35% 


















3.60 


4 


260.17877 


57.4404 


1.10271 


0.21903 




3.60 


5 


259.2686 


•56.79395 


1.104505 


0.218845 




3.60 


' 6 


258.5186 


56.3252 


1.105599 


0.218739 




3.60 


7 


257.9483 


•56.0098 


1.106306 


0.218673 




3.60 


8 


257.6163 


-55.380567 


1.106764 


0.218632 




% Dill. 7/8 




0.13% 


t.i4% 


0.04% 


0.02% 





TABLE 2.1 PHASE CALCULATION VS CYCLE NUMBER 



series of h/2b values and the percent difference for the panel 
code to Theodorsen was plotted in Figure 2.4. This chart shows 
that the h/2b value chosen has a tremendous impact on the code 
results. An h/2b value of .01 gave an acceptable error of 10% 
at Kp = 8. Runs were completed with a value of .01 h/2b and 
the favorable results are shown in Figures 2.53-56. 
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Comparison of the Effect of h/2b Values on Code 
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.2 1 DEGREE PITCH 
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Cl keol 



Real Part of Cl for Panel and Theodorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 
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Kp - 2bw|U 



Figure 2.5 1 Degree pitch Cj^ Re 
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m\ 10 



Figure 



Real Part of Cl for Panel and Theodorsan (pitch, 1.0 dag, .37c, 
NACA0007, 75 panals top and bottom, 3cycles of 65 calc, par cycle) 
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Imag. Part of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 
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Figure 2.7 1 Degree pitch Im 
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Imag. Part of Cl for Panel and Thaordorson (pitch, 1.0 deg, .37c, 
NACA0Q07, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 
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Figure 2.8 1 Degree pitch Cjj Im 
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Cl Mo(j.. 



Mag. of Ci for Panel and Theordorsen (pitch, 1.0 deg. .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 




Figure 2.9 1 Degree pitch Magnitude 
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Mag. of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, Scycles of 65 calc, per cycle) 




Figure 2.10 1 Degree pitch Magnitude 
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Phase of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 




Figure 2.11 1 Degree pitch Phase 
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Phase of Cl for Panel and Theordorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 




Figure 2.12 1 degree pitch Phase 
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Cm Reol 



Cm Real vs Kp for Panel and Theodorsen (pitch, 1.0 deg, ,37c, 
NACA0007, 75 panels top and bottom, Scycles of 65 calc, per cycle) 




Figure 2.13 1 Degree pitch Cjj Re 
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Cm Real 



Cm Real vs Kp for Panel and Theodorsen (pitch, 1.0 deg, .37c, 
NACA0007, 75 panels top and bottom, 3cycles of 65 calc, per cycle) 




Figure 2.14 1 Degree pitch Cj, Re 
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Cm Imag. 



Cm Imag vs Kp for Panel and Theordorsen (pitch, 1 .0 
deg, .37c, NACA0007, 75 panels top and bottom, 
Scycles of 65 calc, per cycle) 

Kp = 2bw/U 




Figure 2.15 1 Degree pitch Cj^ Im 
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Cm Imag. 



Figure 



Cm Imag vs Kp for Panel and Theodorsen (pitch, 1 .0 
deg, .37c, NACA0007, 75 panels top and bottom, 
3cycles of 65 calc, per cycle) 

Kp = 2bw/U 
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Cm MAG. 



Cm Magnitude vs Kp for Panel and Theordorsen 
(pitch, 1 .0 deg, .37c, NACA0007, 75 panels top and 
bottom, Scycles of 65 calc, per cycle) 




Kp = 2bw/U 



Figure 2.17 1 Degree pitch Cjj Magnitude 



40 



Cm MAG. 



Cm Magnitude vs Kp for Panel and Theodorsen (pitch, 
1.0 deg, ,37c, NACA0007, 75 panels top and 
bottom, Scycles of 65 calc, per cycle) 




Figure 2.18 1 Degree pitch Cjj Magnitude 
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Cm Phase 



Cm Phase vs Kp for Panel and Theodorsen (pitch, 1 .0 
deg, .37c, NACA0007, 75 panels top and bottom, 
3cycles of 65 calc, per cycle) 




Figure 2.19 1 Degree pitch Cjj Phase 
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Cm Phase 



Cm Phase vs Kp for Panel and Theodorsen (pitch, 1 .0 
deg, .37c, NACA0007, 75 panels top and bottom, 
Scycles of 65 calc, per cycle) 




Figure 2.2 0 1 Degree pitch Phase 
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TABLE 2.4 6.7 DEGREE PITCH COMPARISON 
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Cl Reol 



Cl Real vs Kp Panel Code (Pitch, 6.7 deg, .37c, NACA 0007, 50 panels 

top and bottom) 




Figure 2.21 6.7 degrees pitch Cj^ Re 
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Cl Real vs Kp Panel Code (Pitch, 6.7 deg, .37, NACA 0007, 75 panels 

top and bottom) 




Kp = 2bw/U 



Figure 2.22 6.7 degrees pitch Re 
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Cl Imag. vs Kp Panel Code (pitch, 6.7 deg, .37c, NACA 0007, 75 panel 

top and bottom) 




Figure 2.23 6.7 degrees pitch Im 
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Cl Imag. vs Kp Panel Code (pitch, 6.7 deg, .37c, NACA 0007, 75 panel 

top and bottom) 




Kp - 2bw/U 



Figure 2.24 6.7 degrees pitch Im 
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Cl Magnitude 



Cl Magnitude vs Kp Panel Code (pitch, 6,7 deg, ,37c, NACA 0007, 50 
panels top and bottom) 




Figure 2.25 6.7 degrees pitch Magnitude 
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Cl Moq 
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Figure 2 



Mag. of CL for Panel and Theodorsen (Pitch,6.7deg, .37 c, NACA 
0007, 75 panels top and bottom, 3cyc65calc.) 
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.26 6.7 degrees pitch Magnitude 
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Cl Phase vs Kp (pitch, 8.7 deg, .37c, NACA 0007, 75 panel top and 

bottom) 




Figure 2.27 6.7 degrees pitch Cj^ phase 
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Cl Phase vs Kp (pitch, 6.7 deg, ,37c, NACA 0007, 50 panel top and 

bottom) 




Kp = 2bw|U 



Figure 2.28 6.7 degrees pitch Phase 
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TABLE 2.5 6.7 DEGREES PITCH COMPARISON 
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Cm Real vs Kp Panel Code (pitch, 6.7deg, ,37c, NACA 0007, 50 
panels top and bottom) 




Figure 2.29 6.7 Degrees pitch Re 
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Cm Reol 



Cm Real vs Kp Panel Code (pitch, 6.7deg, .37c, NACA 0007, 50 
panels top and bottom) 




Figure 2.30 6.7 Degrees pitch Cjj Re 
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Cm Imaginary 
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Figure 2 



Cm imaginary vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 
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Cm Imaginary vs Kp Panel Code (pitch, 6.7deg,.37c, 
NACA 0007, 50 panels top and bottom) 

Kp = 2bw/U 
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Figure 2.32 6.7 Degrees pitch C^j Im 
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Cm Magnitude 



Cm Magnitude vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 




Figure 2.33 6.7 Degrees pitch Cj^ magnitude 
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Cm Magnitudd 



Cm Magnitude vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 




Figure 2.34 6.7 Degrees pitch Cj^ magnitude 
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Cm Phase vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 

0 2 4 6 8 




Figure 2.35 6.7 Degrees pitch phase 
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Cm Phase vs Kp Panel Code (pitch, 6.7deg, .37c, 
NACA 0007, 50 panels top and bottom) 




Figure 2.3 6 6.7 Degrees pitch phase 
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6 PLUNGE h/2b=.0833 COMPARISON 



62 



Real Part of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.37 Plunge h/2b=.0833 Re 
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Real Part of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




0 0.5 1 1.5 2 2.5 

Kp = 2bw/U 



Figure 2.38 Plunge h/2b=.0833 Re 
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Imag Part of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 









— 




*anel 


\ ^ ^ 

\ - 
\ 

\ 

\ 

\ 


Theodorsen 

1 


1 




\ 

\ 

\ 

N 


\ 








\ 



0 2 4 6 8 

Kp = 2bw/U 



Figure 2.39 Plunge h/2b=.0833 Im 

Jj 



65 



CL Imeg 



Imag Part of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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Figure 2.40 Plunge h/2b=.0833 C. Im 

1j 
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Magnitude of Cl for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.41 Plunge h/2b=.0833 Magnitude 
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68 



CL Phaso 



Phase of Cl for Panel and Theodorsen (plunge, .0833 
h/2b, .37c, NACA0007,75 panels top and bottom, 
3cycles of 65 calculations) 




Figure 2.43 Plunge h/2b=.0833 phase 
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CL Phase 



Figure 



Phase of Cl for Panel and Theodorsen (plunge, .0833 
h/2b, .37c, NACA0007.75 panels top and bottom, 
3cycles of 65 calculations) 
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2.44 Plunge h/2b=.0833 a phase 
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Real Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.45 Plunge h/2b=.0833 Re 
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Real Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.46 Plunge h/2b=.0833 C^j Re 
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CM tmag. 



Imag. Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.47 Plunge h/2b=.0833 C„ Im 
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CM Imag. 



Imag. Part of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Kp = 2bw/U 



figure 2.48 Plunge h/2b=.0833 c" 

M 



Im 



CM MAG. 



Magnitude of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.49 Plunge h/2b=.0833 Cjj Magnitude 
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CM MAG. 



Figure 



Magnitude of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 
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CM Phasa 



Phase of CM for Panel and Theodorsen (plunge, 
.0833 h/2b, .37c, NACA0007,75 panels top and 
bottom, 3cycles of 65 calculations) 




Figure 2.51 Pl\mge h/2b=.0833 Cjj Phase 
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8 PLtJNGE h/2b=.01 C„ COMPARISON 
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Real Part of CM for Panel and Theodorsen (plunge, 
.01 h/2b, .37c. NACA0007.100 panels top and 
bottom, Scycles of 65 calculations) 



9 

8 



7 



6 




2 



1 

0 



-1 

0 2 4 6 8 

Kp = 2bw/U 




Figure 2.53 Plunge h/2b=.01 Cjj Re 
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CM Imag. 



Imag Part of CM for Panel and Theodorsen (plunge, .01 
h/2b, .37c, NACA0007,100 panels top and bottom, 
3cycles of 65 calculations) 




Figure 2.54 Plunge h/2b=.01 Cj^ Im 
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CM Phase 



Phase of CM for Panel and Theodorsen (plunge, .01 
h/2b, .37c, NACA0007,100 panels top and bottom, 
Scycles of 65 calculations) 




Figure 2.56 Plunge h/2b=.01 Cj^ Phase 
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Comparison of CM imag vs Kp for Various Values of Pitch 
(NACA0007,100 paneis.l^and 6.7 deg. phch^cyc65)(0c,) 




Figure 2.57 Pitch Oc. 1.0, 3.0, and 6.7 degrees vs Theodorsen 
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1.0, 3.0, and 6.7 degrees 
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III. FLUTTER DETERMINANT 



The proven accuracy of the UPOT Code enabled it to be used 
for the solution of the flutter determinant. 

A. FLUTTER THEORY 

In order to analyze the phenomenon of flutter, it is 
necessary to obtain the equations of motion of the system. To 
simplify the problem the assumption is made that the actual 
motion of the system can be considered a combination of 
fundamental wing bending, and fundamental wing torsion. The 
system can then be replaced by an equivalent system containing 
an airfoil section of unit span restrained by springs against 
independent vertical motion (bending) , and torsion as 
illustrated in Figure 3.1. This paper will not consider the 
aileron hinge case so S and c are set equal to zero. According 
to the class notes of M. Platzer [ref.l] the formulation 
proceeds as follows: 

Consider the balance of the elastic, inertial and 
aerodynamic forces on a mass element: 

• Total Inertial force: -/ dm(h"+rQ'” ) =- (Mh"+S„o!" ) 

Mass : M=/ dm 

Static Moment about the elastic axis: = / rdm 

• The moments about the elastic axis are: 



87 



-!r ih"+ra") dm = - (I^a" +S^h") 



( 3 . 1 ) 



Mass moment of inertia about elastic axis 
= Jr^dm 



• Elastic restoring forces are: -hC^ 

The Equations of motion therefore become: 



h"M+a"S^+hCfj = L 
a"l+h"S+aC„ = M 

a a a 



( 3 . 2 ) 



Where: = Torsional stiffness of the wing 

Cj^ = Stiffness of the wing in translation (plunge) 
M = Mass of the wing per unit span 
These equations can be written in a different way by 
expressing the spring constants in terms of the natural 
frequencies. Consider the airfoil to be so restrained that 
only one degree of freedom is permitted. The equations of 
motion become: 



Mh'^+hC^ = 0 
I„a"l,+aC = 0 

a a a 



so that 



so that 





( 3 . 3 ) 



Hence: = Mcoj,^ C„ = 

The small structural damping of metal aircraft may be 
approximated by a force that opposes the motion and is in 
phase with the velocity. One assumes therefore that the 
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magnitude of the damping is proportional to the elastic 
restoring force. Since the motion of the airfoil is harmonic 
at the critical flutter condition, the structural damping can 
be accounted for by replacing the terms: 
hC,, with hC^(l + ig^) 
ofC^ with q;C|j( 1+ig^) 

Where g^ and g^^ are damping constants multiplied by i to ensure 
that the damping force is in phase with the velocities in the 
simple harmonic motion. 

From equation 3.3 with: 

h(t) = he’“^ and Q!(t) = ae’"^ 

We have : 

h" = (i)^he’“^ and a" = -(o^ae’"^ 

And the equations of motion become: 

=L (3.4) 

( -G)^a J„-<o^hS„+aC„) =M (3.5) 

The equations for the aerodynamic forces were given by 
Fung [ref. 5] and are shown here: 

L=7ipb3(o2 +[!..-( 3 '“'" < 3 . 6 ) 

JD Z* 

Equating equation 3.4 to 3.6 and 3.5 to 3.7 yields: 
Substituting into equation 3.8 and 3.9 for C,^ and and 
using the following dimensional terms: 
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(3.7) 



M=npb‘^u)^ ( [Mf^- ( — +a) L^] — + 



[M- (^+a) (L„+Mu) +(^+a)^L.] a) 



' 2 ' ^ '2 



{-(t)^Mh-u)^aS^+hC)^) =upb^(o^ (Ljj-^ + [L^- { — +a) Lfj] a) (3.8) 

2 



(-(o^a J„-w^/]S„+aC„) =Kp±>^o)2 ( [Af^- ( — +a) L^] — + 
[A/„-(|+a) (L„+W^) + (-|+a)2Lja) 



(3.9) 



p = 



M 



npb‘ 



“ Mb 






Mb'^ 



(3.10) 



The equations simplify and after bringing all terms to the 
left result in: 



A^+Ba=0 

D 

d 4 +Ea = 0 
b 



(3.11) 



This is a homogeneous equation whose solution is obtained 
if the flutter determinant is zero. 



A B 

D E 



= 0 



(3.12) 



Where : 
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(3.13) 



A = )]+L^ 

0) 0)^ n /J 

B=^X„H-L„-L^(y2+a) 

D=jjix„+V 2-L^ (Va+a) 

E=jJir „2 [1- ( (l+igr^) ] -V2 (V2+a) (V2+a) +L^ C/2+a) ^ 

(*) 

^ is the ratio of the mass of the wing to the mass of a 
cylinder of air of a diameter equal to the chord of the wing. 
<i)^ and Wj, are the natural angular frequency (rad/sec) of 
torsional vibration around "a" (elastic axis) and the natural 
frequency in deflection, respectively, x^ is the location of 
center of gravity of the wing measured from a. w is the 
circular frequency of wing vibration. 

The relationships between the code and Theodorsen derived 
earlier in Chapter II can be used here to simplify the 
equations: (note: no damping in this case g„=gh=0) 

For A: manipulating equation 2.20, 






2C 



Lh 






<3.14) 



resulting in: 



0 ). 



= (i[l-(^) 

( 1 ) 



fa)h 

0 ).. 



) ] + 



2C, 



La 






For B: manipulating equation 2.14, 



(3.15) 
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L„-LjV2+a) = 



(3.16) 



KK^a 



resulting in: 



S=^ix„ + 



4C, 



La 



a o 

nK^a 



For D: manipulating equation 2.28 



L„-L^(V 2 +a) =. 



4C, 



Mh 






(3.17) 



(3.18) 



resulting in: 



D=HX„ + 



4C, 



Mh 






(3.19) 



For E: manipulating equation 2.24 



-V 2 (V 2 +a) +M^-L^ C/i+a) +L^(V^+a)^ = 



2 _ 



anKj^ 



(3.20) 



resulting in: 

£;=^ri [l-(.^)^ -t- (3.21) 

w anK/ 

The determinant is expanded to AE-BD=0, and the real and 
imaginary parts are set equal to zero. Substituting {(a}^/o})^= 
X and solving the real (2 roots) and imaginary (1 root) 
equations for values of X corresponding to each reduced 
frequency value. These X values can be plotted as SQRT(X) 
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against Kp and any intersections of real and imaginary parts 
signify a flutter point. 

Knowing that : 



, _ 2b(i) 

U 



and spC=( — ) 

CO 



solve for U , 

critical 



u. 



cri tical 



2b(o^ 

Kpypi 



which is the critical flutter speed. 



(3,22) 



(3.23) 



B . UPOTFLUT CODE 

1. FORMULATION AND INPUT 

The ecjuations derived in the flutter theory section 
above were programmed into a FORTRAN subroutine and attached 
to the UPOT.f code. The UPOT code was modified first to 
enable it to conduct a frecjuency sweep of pitch and plunge 
simultaneously. The resulting frecjuency sweep pitch and 
plunge array data is then sent to the flutter subroutine which 
provides the values of SQRT(x) and Kp for plotting. The 

program also gives a best guess for the U^riticai based on the 
difference betweeen the real and imaginary SQRT (X) values. 
The input file UPOTFLUT. IN is very similar to the regular 
UPOT. IN file with the addition of actual physical properties 
of the system being analyzed. The user should start the 



93 



analysis in the pitch mode first (lOSCIL =1, ITRANS =0) to 
ensure complete coverage of all frequencies of interest. The 
following relations were taken from NACA TR-685 [ref. 8] and 
should prove helpful in determining the physical properties 
needed for program operation. 

K = mass ratio = 7rpb^/M 
K = 1//Z 
X„ = S„/Mb 
r„2 = i^/Mb2 
2 . OUTPUT 

Outputs from the code have been limited to reduce the 
amount of computer space taken up by the code operation. A 
sample input and output file are contained in Figure 3.2. 
The following list describes the input/output files and the 



data 


they contain: 






a. 


UPOTFLUT.IN 


The 


input file 


figure 3.2a 


b. 


CL.d 


Same 


as UPOT.f 


output 


c . 


CM.d 


Same 


as UPOT.f 


output 


d. 


PHASE, d 


Same 


as UPOT.f 


output 


e . 


CPSS.d 


Same 


as UPOT.f 


output 


f . 


CPU005.d 


Same 


as UPOT.f 


output 


g- 


PHZSWP . d 


This 


file contains 


h. 


PLHZSWP . d 


This 


file contains 


i . 


PITCH. in 


This 


file contains Kp , Re, 

Im 
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b = semichord (ft.) 

cb = distance between midchord and aileron hinge, positive if aft of 
midchord (ft.) 

eb *= distance between midchord and aileron leading edge, positive 
■* aft of midchord (ft.) 

= distance between rotation point (elastic axis) and midchord, 
positive if aft of midchord (ft.) 

h = bending deflection of rotation point (elastic axis), positive down- 
ward (ft.) 

, a = angular deflection about rotation point (elastic axis), positive 
for leading edge up (radians) 

P == angular deflection of aileron about aileron hinge relative to wing 
chord, positive for aileron leading edge up (radians) 



Figure 3.1 Simplified System Geometry 
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$tdln 



Page 1 



AIRFOIL TYPE ! NACA 0012 AIRFOIL 

NLOWER = 50 , NUPPER = 50 



IFI.AG 


NLOWER 


NUPPER 




0 


50 


50 




AIRFOIL TYPE 
1 






IRAMP 


lOSCIL 


ALPI ALPMAX 


PIVOT 


0 


1 


-3.0 3.0 


0.3 


FREO 


RFOSTP 


rfofnl 




.85 


.01 


,95 




IGUST 


UGUST 


VGUST 




0 


0. 


0. 




ITRANS 


DELMX 


DELHY DELI PHASE 




0 


0.00 


.0833 -.0833 0.00 




CYCLE 


NTCYCLF 


TOI. 




3 


65 


0.005 




naot i 


naot X 


aoa values multiplied 


by 10 { 


2 


05 10 


20 25 39 50 





Semi-chord Walpha Wh Mass 

6 90 72.0 .537R9 

lalpha Salpha Density 

4.84102 .645468 .002378 

Comments . . . 



TRAMP 0: n/a 



RFREO Is based on full chord 



1: Straight ramp 
2: Modified ramp 



lOSCIL 0: n/a RFREO is based on full chord 

1: Sinusoidal pitch, motion starts at min Aoa 

ITRANS 0: n/a 

1: Trans 1 at lonal harmonic oscillation 

CYCLE : I of cycles for oscillatory motions 

-In case of ramp, cycle=l,5 denotes airfoil Is held 
at max aoa for the duration of .5 cycle 
-For steady state solution set It to 0 

NTCYCLE; f of time steps for each cycle 

CYCLE*NTCYCLE Is limited to 200 currently. 

NAOT : I of input aoa for cp output 

- angles should be in increasing order, 

- for oscilatory motions angles should increase 
first, then decrease. Decreasing angles are for 
the return cycle., 

SEMI-c5hORD Half Chord in feet. 

Walpha,Wh, uncoupled natural frequencies of the system in question. 
Walpha is pitch and Wh is piunge(HZ). 

Mass specific mass of the system In slugs/foot of span 

lalpha Moment of Inertia of system about the elastic axis(a) 

per unit span length. 

Salpha Static moment of wing-aileron per unit span length 

Density Mass of air per unit of volumefslugs per ft''3) 



Figure 3.2a UPOTFLUT.in example input file 
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ijtdln 


Page ^ 




AIRFOIL TYPE : NACA 0012 AIRFOIL 
NLOWER = 50 , NUPPER > 50 


OSCILLATORY MOTION, lOSCIL 


1 




FREO SWEEP 
FREO = 0.850000 








PHASE SHIFT ANALYSIS 






FREO = 


0.8500000 






w 0.8500000 

kp^ 0.8500000 Ifreq 


1 






AMPLITUDE? clamp, cmamp 


0.2316691 


2 .9221054E-02 




ioscll » Oitrans = 0 






PHASE; clp, cmp 


202.4033 


-63.70797 




AVERAGE DRAG, TOTAL DRAG 


1.5930884E-03 


0. 1051438 




ETAS, WBAR 


-0.1430808 


-1 . 11 34184E-02 




PHASE SHIFT ANALYSIS 






FREO = 


0.8500000 






w 0.8500000 

kp= 0.8500000 Ifreq 


1 






AMPLITUDE; clamp, cmamp : 


0.2779089 


2.2117507E-02 




ioscll = Oitrans = 1 






PHASE; clp, cmp 


270. 5030 


37.10253 




AVERAGE DRAG, TOTAL DRAG 


-6.2051453E-03 


-0.4095396 




ETAS, WBAR 


0.1050769 


-0.1181067 




FREO SWEEP 
FREO = 0.860000 








PHASE SHIFT ANALYSIS 






FREO = 


0.8600000 






w 0.8600000 

kp= 0.8600000 Ifreq 


2 






AMPLITUDE? clamp, cmamp ; 


: 0.2317267 


2.9490557E-02 




ioscll = Oitrans = 0 






PHASE; clpr c:mp 


202.8760 


-63.85251 




AVERAGE DRAG, TOTAL DRAG : 


: 1.5921656E-03 


0.1050829 




ETAS, WBAR ! 


r -0.1398697 


-1.1383206E-02 




PHASE SHIFT ANALYSIS 






FREO = 


0.8600000 






w 0.6600000 

kp= 0.8600000 ifreq 


2 






AMPLITUDE? clamp, cmamp ; 


r 0.2804866 


2.2494521E-02 




ioscil = Oitrans = 1 






PHASE; clp, cmp 


270.7843 


36.79393 




AVERAGE DRAG, TOTAL DRAG 


-6.3188463E-03 


-0.4170439 




ET?^S, WBAR 


0.1047920 


-0.1205979 




FREO SWEEP 
FREO = 0.870000 








PHASE SHIFT ANALYSIS 






FREO = 


0.8700000 






w 0.8700000 









Figure 3.2b UPOTFLUT exainple output file 
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X ^ V 


■ ^ -Page 2^ 


kp= 0.8700000 Ifreq 


3 


AMPLITUDE; clamp, cmamp 


0.2317936 2.9765518E-02 


loscil = Oitrans 0 


PHASE; clp, cmp 


203.3604 -64.02047 


AVERAGE DRAG, TOTAL DRAG 


1.5893428E-03 0.1048966 


ETAS, WBAR 


-0.1365460 -1 .1639614E-02 


PHASE SHIFT ANALYSIS 


FREO = 


0.8700000 


w 0.8700000 




kp= 0.8700000 Ifreq 


3 


AMPLITUDE; clamp, cmamp 


0.2830637 2 . 2874046E-02 


loscil - Oitrans = 1 


PHASE; clp, cmp 


271.0655 36.47753 


AVERAGE DRAG, TOTAL DRAG 


-6.4334869E-03 -0.4246101 


ETAS, WBAR 


0.1045149 -0.1231114 


FREO SWEEP 




FREO = 0.880000 




PHASE SHIFT ANALYSIS 


FREO = 


0.8800000 


w 0.8000000 




kp- 0.0000000 Ifreq 


4 


AMPLITUDE; clamp, cmamp ! 


: 0.2318744 3.0037180E-02 


loscil = Oitrans = 0 


PHASE; clp, cmp 


: 203.8330 -64.15525 


AVERAGE DRAG, TOTAL DRAG : 


; 1 .5880425E-03 0.1048108 


ETAS, WBAR ! 


r -0.1335113 -1 . 1894441E-02 


PHASE SHIFT ANALYSIS 


FREO = 


0.8800000 


w 0.8800000 




kp= 0.8800000 Ifreq 


4 


AMPLITUDE; clamp, cmamp i 


t 0.2856959 2 . 32 582 35E-02 


loscil - Oitrans = 1 


PHASE; clp, cmp 


271.3546 36.17284 


AVERAGE DRAG, TOTAL DRAG 


-6.5488103E-03 -0.4322215 


ETAS, WBAR 


0.1042214 -0.1256711 


FREO SWEEP 




FREQ = 0.890000 




PHASE SHIFT ANALYSIS 


FREQ = 


0.8900000 


w 0.8900000 




kp=^ 0.8900000 ifreq 


5 


AMPLITUDE; clamp, cmamp ; 


0.2319494 3.0311935E-02 


loscil = Oitrans = 0 


PHASE; clp, cmp 


204.3115 -64.29196 


AVERAGE DRAG, TOTAL DRAG 


1.5859993E-03 0.1046759 


ETAS, WBAR 


-0.1304959 -1 .2153637E-02 


PHASE SHIFT ANALYSIS 


' FREO = 


0.8900000 


w 0.8900000 




kp= 0,8900000 Ifreq 


5 


AMPLITUDE; clamp, cmamp : 


: 0.2883391 2 . 3 64 4 408E-02 



Figure 3.2c UPOTFLUT example output file 
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■ stdln ■ piggi": 


loscll = Oltrans = 1 


PHASE; clp, cmp 


271.6436 35.87401 


AVERAGE DRAG, TOTAL DRAG 


-6.6649993E-03 -0.4398900 


ETAS, WBAR 


0.1039310 -0.1282581 


FREQ SWEEP 




FREO = 0.900000 




PHASE SHIFT ANALYSIS 


FREQ = 


0.9000000 


w 0.9000000 




kp= 0.9000000 Ifreq 


6 


AMPLITUDE; clamp, cmamp 


0.7320280 3.0592278E-02 


loscll = Oltrans =• 0 


PHASE; clp, cmp 


204.7803 -64.43259 


AVERAGE DRAG, TOTAL DRAG 


1 .5828811E-03 0.1044701 


ETAS, W8AR 


-0.127461 7 -1 .241 8402F.-O7 


PHASE SHIFT ANALYSIS 


FRFQ ^ 


0.9000000 


w 0.9000000 




kp» 0.9000000 Ifreq 


6 


AMPLITUDE; clamp, cmamp 


0.2909827 2 . 4039967E-02 


loscll = Oltrans = 1 


PHASE; clp, cmp 


271.9288 35.58691 


AVERAGE DRAG, TOTAL DRAG 


-6.7821071E-03 -0.4476191 


ETAS, WBAR 


0.1036481 -0.1308680 


FREQ SWEEP 




FREQ = 0.910000 




PHASE SHIFT ANALYSIS 


FREQ = 


0.9100000 


w 0.9100000 




kp= 0.9100000 ifreq 


7 


AMPLITUDE; clamp, cmamp ; 


: 0.2321466 3 . 08651 42E-02 


loscll = Oltrans = 0 


PHASE; clp, cmp 


205.2529 -64.61227 


AVERAGE DRAG, TOTAL DRAG 


1.5800473E-03 0.1042831 


ETAS, WBAR 


-0.1245367 -1 .2687406E-02 


PHASE SHIFT ANALYSIS 


FREQ = 


0.9100000 


w 0.9100000 




kp= 0.9100000 Ifreq 


7 


AMPLITUDE; clamp, cmamp 


: 0.2936268 2 . 4 4404 34E-02 


loscll = Oltrans = 1 


PHASE; clp, cmp 


272.2159 35.29980 


AVERAGE DRAG, TOTAL DRAG 


-6. 9000078E-03 -0.4554005 


ETAS, WBAR 


0.1033707 -0.1335003 


FREQ SWEEP 




FREQ = 0.920000 




PHASE SHIFT ANALYSIS 


‘ FREQ = 


0. 9200000 


w 0.9200000 




kp= 0.9200000 ifreq 


8 



Figure 3. 2d UPOTPLUT example output file 
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AMPLUTJDf ; cl imp, cmawp : 0.7177 7S9 1 . 1 nft^97r-0? 

loac II* OltfADj- 0 

PHASC; dp, rmp : 205.7313 -t4.«<»62« 

AVERAfTE DRAO, TOTM ORAr; : 1 . 5 782 3 72F-03 0.104U37 



r.TA5:, WflAR 



0.9700000 
Vp- 0.9700000 



>0. 1718783 



PHASE Siliri ANALYSIS 
FPro - 0.9700000 



I f 



AMPLITt’Or.; rl.imp, cmainp ; 0.7987 705 

lotcll = Oltr.ins * 



-1 .2949497F.-0? 



?.4fl435e7F-07 



PHASE; cip, rinp 

AVERAGE DRAG, TOTAI DRAG 
FTAS, wnAR 



rPEO swrrp 
FRro - 0.930000 



2 72.501 1 



35.01288 



-7.Ol0803Or.-O3 -0.4832331 
0.1030«B7 -0.1381548 



PHASE SHiri ANALYSIS 
7 RfO - 0. 9799999 



Ifr-T 



w 0.9799999 
kp- 0.9799999 



AMPLlUlOr,; clump, cmamp t 0.7374548 
loscl I’' Oltranj- 0 

PHASE; «'lp» ^»"P ; 208.2000 



3. 14173IOr-07 
-84.02713 



AVERAGE DBAG, TOTAL DRAG t 
F.TAS, WRAR t 



1.5755824r-03 0.1039071 



-0.1191857 



PHASE SHTPT ANALYSIS 
f RPO • 0. 9799999 



1 f Tcq 



w 0.9799999 
kp- 0.9799999 



AMPLITDDP; clamp, cm«mp i 0.7989171 
loscl] - Oltrans * 



-1 . 3721 807P.-07 



7.5750337P-07 



PHASE; dp, cmp 

AVERAGE DRAG, TOTA7, DRAG 
F.TAS, WBAB 



FREO SWEFP 

FRF.O • 0.940000 



2 72.7 901 



34.72940 



-7.1302 1O0F.-O3 -0.4711219 
0.1078378 -0.1388317 



w 0.9399999 
kp- 0.9399999 



FHASE SHIFT ANALYSIS 
FBFO - 0.9399999 



1 f rcq 



10 



AIIPLIIUDF; cl. amp, cmsmp : 0.7 3784 74 

losdl Oltrans • 

PHASE; dp, cmp : 208.8048 

AVERAGE drag. TOTAL DRAG : 1 , 5 7 15322F-03 0.1037211 

F.TAS, WflAR : -0.1 1 84358 -I . 34989R9E-0? 



1.1 8990 78F.-0? 

-84.94021 



PHASE SH17T ANALYSIS 

FRro - 0, 9199999 



! f rcq 



0. 9399999 
Vp-T 0.9304999 



AMFLITODF.; clamp, cmamp : 0.3035810 

loac 1 1 • Olt ran * ^ 



7 .5859987F-07 



PHASE; <*lp» ^■mp 

average drag, TO! At. DRAG 
FIAS, WflAR 
Humber of Kp Values 
Plunge Valuc/Fiill CherH Oi/7hl 
Alpha 

PIVOT POINT (a) Of Flaatlr Axis 
Half Chord (b) 

Alpha Uncoupled Nat. 7req. 
plunge Oncouple/4 War. Freq 
1 alpha 
S alpha 
has a 

Alt Density 
Hass Ratio 
Dist. To Wing CG Aft of Elas. Ax1s> 



34.44823 



273.0792 
-7.2507137F-O3 -0.4790751 
0.1075748 -0.1415304 

10 

0.0033 
3.0000 
-0.4000 
€.0000 
♦ 0.0000 
77.0000 

4.8410 
0.8455 
0.5379 
0.0074 

2.0000 

0.2000 



Rp crll a 0.9000000 
Dirr * 1 . 1507273f-03 

SORTX » 0.9579420 

U cr It - I7S7 . 881 



Figure 3.2e UPOTFLUT example output file 
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j. PLUNGE. in This file contains K_, RE, Cj^ RE, 

CMh IM 

k. FLUTPLOT.d This file contains K^,SQRT(x) Re, SQRT(X) RE, 

SQRT(X) IM 

3 . VALIDATION 

The program was tested against some sample cases to 
check for code validity. The first case was taken from 
reference 6 example #1, p. 236. Figures 3.3 and 3.4 show 
plots of the FLUTPOT.d file. Figure 3.3 shows the initial 
look over a wide range of Kp and after finding the approximate 
flutter location Figure 3.4 shows a closer look at the Kp 
range of interest. This example calculated a Ucritical 
161.985 ft/sec. which compares favorably to the example value 
of 162 ft/sec. The next example was taken from NACA TR-685 
[ref. 8] case #1 p. 8. Figures 3.5 and 3.6 again show the 
initial and final looks for this analysis. The example called 
for a U^j-itical miles/hr and the program returned a 

value of 570 miles/hr. Next, the code was tested over a range 
of ratios as done in NACA TR-685, pll., graph I-A(a) . 

Figure 3.7 shows the comparison between the two methods. 
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Figure 3.3 Initial look at flutter, excuonple 1 
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Figure 3.4 Final look at flutter, example 1 
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Figure 3.5 Initial look at flutter, example 2 
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Figure 3.6 Final look at flutter, example 2 
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Figure 3.7 Calculations 
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IV. FLOW VISUALIZATION EXPERIMENT 



A. INTRODUCTION 

The purpose of this experiment was to document the 
production of thrust by a plunging airfoil. This was a 
preliminary experiment to better understand the vortex pattern 
produced by a plunging airfoil, and to examine the production 
of thrust using smoke flow visualization techniques. 

An explanation of what constitutes a propulsive vortical 
signature along with smoke flow visualization of the 
propulsive vortical patterns is given in Reference 7. In this 
reference, the explanation is given by contrasting the 
vortical pattern produced by a cylinder (drag) with the 
vortical pattern produced by a plunging airfoil (thrust) . The 
cylinder produced a vortical sheet where the top row of 
vortices rotated clockwise and the bottom row of vortices 
rotated counterclockwise. This pattern induces a velocity 
component in the upstream direction (Biot-Savart law) . In 
contrast, the plunging airfoil produced a clockwise rotating 
vortex sheet on the bottom row. This pattern induces a 
velocity component in the downstream direction. Reproduction 
of the flow visualization data from Reference 7 is shown in 
Figure 4.1. 
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Figure 4.1 View of flow over cylinder (top) and plunging 
airfoil (bottom) [Ref. 7] 
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B . THEORY 



A comparison was done using the incompressible panel code, 
U2DIIF. The purpose of this study was to examine the vortical 
pattern produced by the panel code, and determine if the 
vortical signature matched experimental results. The input to 
the panel code was set up to best match the conditions of the 
experiment described in the next section. The panel code was 
run using a plunge amplitude, h/2b equal to .1977, a reduced 
frequency of 1.8 and a zero mean AOA. The results of the 
vortical pattern are shown in Figure 4.2. Aside from the 
starting vortex, this is clearly a thrust producing vortical 
sheet. Furthermore, the vortical pattern is similar to that 
produced by the experiment shown in Figures 4.10 and 4.11. 

C . EXPERIMENTAL SETUP 
1. Plunging Airfoil 

The plunging airfoil used in this paper was originally 
a wing taken from the rotor of a model helicopter. The wing 
was attached to a MB250 Shaker Table as shown in Figure 4.3 
The wing was made from a NACA0007 airfoil section and 
consisted of a 2.45" chord and a 22" span. The wing was built 
from a foam core and finished with a layer of graphite epoxy 
composite for added fatigue strength. The airfoil's drive 
mechanism was a MB 250 Shaker Table capable of 1" total 
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Figure 4.2 Wake pattern produced by U2DIIP code 
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Figure 4.3 Shaker table with wing 
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deflection. The Shaker Table was limited only by resonance 
frequencies of the wing which occurred around 20 Hz or 1200 
rpm. 

2. WIND TUNNEL 

The wind tunnel used in this experiment was a very low 
speed, low turbulence smoke tunnel. It is made of plexiglass 
walls and a contraction ratio of 2.8:1. The motor provides 
wind tunnel velocities between 0 and 10 feet per second (fps) . 
The smoke was created using a Rosco smoke generator and piped 
into the tunnel in the test section using a small seven tube 
smoke rake constructed for this experiment. Figure 4.4 is a 
photograph of the wind tunnel and smoke rake used in this 
experiment . 

D. TEST PROCEDURE 

Testing was conducted in the low speed smoke tunnel under 
several different conditions. The speeds of the tunnel were 
approximately 1.04 fps, 1.47 fps, and 1.56 fps (measured 
visually) . These low speeds allowed good pictures and the 
ability to get higher reduced frequencies without calling for 
too high a load on the wing. The actual plunging harmonic 
frequencies ranged from 1 to 15 Hz and amplitudes from 1/16" 
to 1" peak to peak. The tunnel was initially turned off and 
the Shaker Table turned on with stagnant smoke in the tunnel. 
The purpose was to see if the 

plunging airfoil would draw the smoke through the tunnel like 
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Figure 4.4 Tunnel and smoke rake 
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a fan, thus showing the production of thrust by the plunging 
airfoil . 

Photos were taken using a Nikon 3 5inm camera and Kodak 
TMAX-400 ASA black and white film. The shutter speed was set 
to 1/125 seconds with an aperture setting of 4.0 for the light 
conditions. Film developing time was optimized at 9 minutes at 
75 degrees F. 

E. RESULTS AND DISCUSSION 

The result for the tunnel off condition flow visualization 
experiment was as expected. The wing in fact accelerated the 
smoke in its vicinity. 

The result of the additional rake flow visualization 
experiments are shown in Figures 4.5-4.14. Figure 4.5 shows 
the stationary airfoil at zero degree AOA. The Reynolds number 
(based on airfoil chord) is 10,000. It can be seen that the 
airfoil produces a small wake with the boundary layer mostly 
attached. Figures 4.6 through 4.14 show the vortical wake 
flow patterns produced by plunge oscillations at various 
frequencies as indicated. Most of these pictures reveal the 
propulsive vortical street pattern discovered in Reference 7 . 
Previous experiments by Neace, [ref. 9] found that the tunnel 
was too small for the airfoil size used, but the airfoil size 
for the present experiment seemed to be optimum, as seen by 
the long trail of vortices. The vortical patterns show that 
the bottom vortex is rotating clockwise, and the top vortex is 
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Figure 4.5 Steady airfoil 1.56 ft/s 
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Figure 4.6 




Kp = 1.8008, h/2b=.1977, 1.56 ft/s 









j;i>- 
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Figure 4.7 = 2.467, h/2b=. 10204, 1.56 ft/s 
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Figure 4.8 Kp = 2.467, h/2b=.1913, 1.56 ft/s 
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Figure 4.9 Kp = 4.112, h/2b=. 14031, 1.56 ft/s 
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Figure 4.10 Kp = 6.167, h/2b=. 05102, 1.56 ft/s 
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Figure 4.11 Kp = 6.753, h/2b=.1084, 1.56 ft/s 
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Figure 4.13 =8.223, h/2b=. 01913, 1.56 ft/s 
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Figure 4.14, Kp=12.335, h/2b=. 01775, 1.56 ft/s 
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rotating counterclockwise, which is a thrust producing 
vortical sheet. It can be seen in the pictures that the 
frequency greatly affects the vortical strength (size) . 
Increasing the frequency leads to an increase in wake 
vorticity . 
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V. LIFT ENHANCEMENT PRODUCED BY A PLUNGING AIRFOIL 



A. THEORY 

Chapter IV demonstrated the propulsive capability of a 
plunging airfoil. The production of thrust implies the 
generation of a jet flow which, in turn, may be utilized as a 
boundary layer control device. Therefore, an additional test 
was conducted in the NPS smoke tunnel in order to explore the 
feasibility of this concept. 

B . SETUP 

The same NACA0007 plunging wing was used as in chapter IV 
with a different driving mechanism. The wing was mounted to 
an ELECTRO- SEIS Model 113 Shaker Table by APS Dynamics, Inc. 
The shaker was located below the test section of the NPS smoke 
tunnel (Figure 5.1). The plunging airfoil was mounted to 
struts at both ends to prevent excessive bending while 
plunging. The large airfoil is a cambered profile taken from 
the rotor of a full size helicopter (13" chord, 2" Thickness, 
modified NACA airfoil, Reynolds Number of 52,000) . The wing is 
suspended from the tunnel ceiling as shown in the flow 
pictures. The design allowed for full movement of the big wing 
to position it in the vicinity of the plunging airfoil. 
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C . WIND TDNNEL 



This study used the Naval Postgraduate School's flow 
visualization wind tunnel. The tunnel is an open- circuit one, 
with air entering an inlet that measures 4.5 m X 4.5 m 
(15'X15') . As the air enters the tunnel, it passes through a 
7. 5 -cm long honeycomb. A 9:1 ratio square contraction cone 
directs the flow into a test section that is 1.5 m X 1.5 m 
(5'X5'), and 6.7 m (22') long, as seen in Figure 5.2. The 
flow is then exhausted into the atmosphere through a fan, 
which uses variable pitch blades to control the speed of the 
flow. The speed control toggle switch is located right below 
the red and green on/off switch located in the left side of 
the tunnel control room. The tunnel speed was determined using 
a digital manometer which was verified for accuracy (Figure 
5.3) . 

An observation booth is located on the side of the tunnel. 
A glass window, 1.6 m X 1.1 m (5.2' X 3.4'), provides the 
primary viewing area from the observation room and a second 
one, 0.4 m X 1.23 m (1.33'X4'), is located in the tunnel's 
roof. The main viewing window had sufficient area for most of 
the photography, with the top window used for illiimination. A 
circular turntable was located on the floor of the test 
section [ref. 11] which allowed for easy access to the shaker 
table. The walls and floor of the test section were flat black 
for low light reflectivity. 
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Figure 5.1 Shaker table setup below tunnel 



28 







TLAF Tliy 




Figure 5.2 Tunnel Layout 
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D. SMOKE GENERATION 



The smoke was generated in the Rosco smoke/fog machine. 
Many different smoke injection techniques were tried but with 
less than satisfactory results. Smoke rakes were first tried 
outside the tunnel with the tube number varying from 2 to 30 
tubes . The tubes were inserted in the honeycomb and also 
separated different distances from the inlet of the tunnel. 
The rake was also tried inside the test section with very bad 
results (smoke dispersed immediately) . Problems ranged from 
lack of smoke and turbulence when enough smoke was present . 
The Rosco machine at its lowest setting was producing a very 
high smoke volume and whenever the flow was restricted by a 
smoke rake the smoke production went way down. The final 
technique adopted was very simple. The smoke output was sent 
directly from the machine to a 1" nozzle which was manually 
waved at the entrance of the tunnel to make a steady cloud. 
The steady cloud was gradually pulled into the tunnel, 
producing a thick smoke sheet in the test section. 

E . PHOTOGRAPHY 

Photos were taken using a Nikon 35mm camera and Kodak 
TMAX-400 ASA black and white film. The film speed was set to 
1/250 seconds with an aperture setting of 4.0 for the light 
conditions. Film developing time was optimized at 9 minutes at 
75 degrees F. 
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F . EXPERIMENTAL PROCEDURES 

The first step involved a look at the large airfoil to 
verify normal flow patterns (Figure 5.4 and 5.5) and find the 
AOA for initial trailing edge separation (Figure 5.6). Next, 
the plunging airfoil was placed in the tunnel by itself and a 
run was made to verify the propulsive capability in the larger 
tunnel at higher speeds. As seen in Figure 5.7, the airfoil 
produced a drag vortical flow. Figure 5.8 and 5.9 shows the 
propulsive pattern of the propulsive airfoil. Finally, the two 
airfoils were placed in close relative position to see the 
interference effect between the two airfoils. 

G. RESULTS AND DISCUSSION 

Several airfoil position combinations were studied, as 
shown in Figures 5.10 through 5.15. Figure 5.10 and 5.11 show 
the plunging airfoil located at approximately .65 chord of the 
large airfoil. Figure 5.12 through 5.15 show the plunging 
airfoil located at approximately .75 chord. 

The differences between the plunging on and off condition 
were not easy to see with the eye but pictures indeed showed 
some differences between the two conditions. A shortcoming of 
this experiment was the inability of the plunging airfoil to 
run parallel with the large airfoil. Additionally, sizing of 
and relative positioning of the two airfoils was not optimized 
to give best results. The two airfoils were chosen from the 
resources available and time constraints prevented a more 
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Figure 5.4 Large airfoil at zero AOA 
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Figure 5.5 Large airfoil at 10 degrees AOA 
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Figure 5.6 Large airfoil at 12 degrees AOA (stall) 
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Figure 5.7 Plunge airfoil steady 
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Figure 5 




8 Plunging airfoil Kp=1.71 
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Figuirs 5.9 Plunging airfoil Kp— 3.42 
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Figure 5.10 Large airfoil at 12 degrees AOA, steady plunging 
airfoil position 1 
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Figure 5.11 Large airfoil at 12 degrees AOA, plunging 
airfoil Kp=3.42, position 1 
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Figure 5.12 Large airfoil at 12 degrees AOA, plunging 
airfoil steady, position 2 
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Figure 5.13 Large airfoil at 12 degrees AOA, plunging 
airfoil Kp=3.42, position 2 
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Figure 5.14 Large airfoil at 14 degrees AOA, plunging 
airfoil steady, position 2 
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Figure 5.15 Large airfoil at 14 degrees AOA, plunging 
airfoil Kp=3.42, position 2 
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detailed investigation of the interference effects between the 
two airfoils. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. SINGLE AIRFOIL ANALYSIS 

The modified version of U2DIIF (UPOT) can perform 
aerodynamic calculations over any range of reduced 
frequencies. The nonlinear theory presented here for harmonic 
motion, and the phase relationships that exist between the 
airfoil motion and the aerodynamic forces have been 
extensively verified by comparison with Theodorsen's linear 
theoiry. Furthermore, this panel code was applied to the 
analysis of incompressible bending- torsion airfoil flutter. 
Again, excellent agreement with the classical Theodorsen 
analysis was obtained. 

Access to faster computational means is recommended to 
shorten the time needed to predict the flutter points. The 
code should be modified to incorporate three-dimensional 
calculations which would help solve more difficult flutter 
problems . 

B. FLOW VISUALIZATION EXPERIMENTS 

The flow visualization experiment successfully showed the 
development of thrust produced by a plunging airfoil. The 
enhanced lift experiment, on the other hand, was not a 
complete success. The smoke visualization presented 
difficulties that were not satisfactorily overcome. As a 
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result, the pictures taken were somewhat inconclusive. 
Futhermore, the angle of attack of the oscillating airfoil 
could not be changed thus making it difficult to achieve a 
flow condition conducive to lift enhancement. 

It is recommended that further experiments be conducted in 
the low speed smoke tunnel with a shaker table capable of 
moving an airfoil at harmonic frequencies near 40 HZ. 
Additionally, the airfoil must be modified to allow change of 
AOA. Finally, the Rosco smoke machine output volume must be 
modified to permit much lower smoke output. This final point 
proved to be the single largest detriment to the visualization 
experiment . 
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